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Part  I 


THE  EXTINCTION  BY  SMALL  ALUMINUM  PARTICLES 


Written  by  Janice  Rathmann 


EXTINCTION  MEASUREMENTS  ON  A!  UMINUM  AND  CARBON  SMOKE  PARTICLES 
FROM  FAR  INFRARED  TO  FAR  ULTRAVIOLET 


1.  ABSTRACT 

Absorption  measurements  were  performed  from  .12um  to 
250ym  on  small  particles  cf  Al  produced  by  the  standard 
inert  gas  evaporation  technique.  These  were  done  because 
experimental  results  of  Granqvist  et.al.  revealed  a 
discrepancy  of  over  three  orders  of  magnitude  with  the 
predictions  of  both  the  classical  Drude  theory  for  spherical 
particles  and  the  quantum  mechanical  Gor'kov  Eliashberg 
theory.  Our  measurements  serve  to  confirm  the  results  of 
Granqvist  et.al.  in  the  far  1R  and  to  extend  the  data  to 
the  UV.  In  the  present  paper  the  "anomalous"  absorption  is 
explained  with  a  model  using  the  classical  Drude  theory  by 
taking  into  account  a  distribution  of  ellipsoidal  shapes  and 
an  amorphous  oxide  layer,  which  is  found  to  cover  the 
particles.  In  order  to  determine  the  optical  constants  of 
amorphous  aluminum  oxide  required  in  the  calculation,  a  new 
technique  applicable  to  powdered  materials  was  used.  Our 
model  with  the  shape  distribution  was  found  to  give 
excellent  agreement  with  the  experimental  results.  Several 
models  proposed  by  other  authors  to  explain  the  "anomalous" 
far  IR  absorption  were  also  considered,  but  did  not  agree  as 


2 .  INTRODUCTION 


The  absorption  of  light  by  several  sizes  of  small 
aluminum  particles  was  measured  in  the  far  infrared  (from  3 
cm"1  to  150  cm-1)  by  Tanner,  Sievers,  and  Buhrman  (1975). 
For  the  smallest  particle  size  they  observed  a  multiple  peak 
structure  which  they  interpreted  in  terms  of  quantum  size 
effects  (QSE)  predicted  by  the  Gor ' kov  Eliashberg  theory 
(1965).  Granqvist  et.al.  (1976)  repeated  the  experiment 
using  an  improved  sample  fabrication  technique  that  produced 
smaller  particles  of  a  more  uniform  size.  The  structure 
that  was  expected  to  occur  due  to  QS.E  was  not  seen  and  the 
absorption  they  measured  was  three  orders  of  magnitude 
higher  than  either  the  Gor ' kov  Eliashberg  theory  or  the 
classical  Drude  theory  predict  (see  Fig.  1.  ) .  The  failure 
to  see  structure  caused  by  QSE  is  understandable  since  the 
size  distribution  of  particles  was  not  sharp  enough 
(Granqvist  et.  al.  1976,  Granqvist  1978).  In  fact  Devaty 
and  Sievers  (1980)  point  out  that  it  is  very  unlikely  that 
QSE  can  be  seen  with  present  particle  preparation 
techniques.  The  anomalously  high  absorption  is  not  as 


LOG  FREQUENCY  (cm'1) 


Absorption  Coefficient  of  Al 

The  data  of  Granqvist  £_fc  al .  (solid 
circles)  and  calculated  absorption 
for  spheres  (dashed  line) 
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easily  explained.  There  have  been  many  theories  proposed  to 
account  for  this,  including  the  following. 

1)  The  aluminum  particles  produced  by  Granqvist 
et.al.  have  an  amorphous  oxide  coating  on  them  (they 
purposely  let  in  oxygen  in  the  inert  gas  evaporation  process 
to  ensure  electrical  neutrality  of  the  particles) .  It  has 
been  suggested  by  Simanek  (1977)  that  the  absorption  occurs 
primarily  in  the  oxide  coating.  His  model  requires  the 
particles  to  aggregate  into  clusters  so  that  in  effect  the 
sample  consists  of  randomly  oriented  cylinders  of  aluminum 
oxide  with  aluminum  particles  embedded  within. 

2)  Glick  and  Yorke  (1978)  suggest  that  a  large  part 
of  the  absorption  in  the  far  infrared  should  be  due  to 
phonon  excitation.  The  Gor ' kov  Eliashberg  theory  and  most 
of  the  others  assume  purely  electronic  absorption.  They  say 
the  structure  of  the  lattice  absorption  should  closely 
parallel  the  structure  in  the  phonon  density  of  states  of 
crystalline  aluminum  and  expect  the  electronic  absorption  to 
have  no  appreciable  structure.  The  absorption  they 
calculate  is  in  better  agreement  with  the  data  of  Granqvist 
et.al.  (as  far  as  it  extends)  than  either  the  Gor' kov 
Eliashberg  theory  or  the  Drude  theory. 

3)  Another  theory  (Lushnikov,  Maksimenko,  and 
Simonov  1978,  and  Maksimenko,  Simonov  and  Lushnikov  1977) 
says  that  when  the  particle  size  is  close  to  or  smaller  than 
the  mean  free  path  of  the  material  there  are  QSE  due  to 


correlations  between  the  one  electron  energy  levels  inside 
the  particle  (Coulomb  interactions  between  conduction 
electrons) .  This  theory  produces  some  improvement  but  does 
not  completely  account  for  tfc  three  orders  of  magnitude 
difference  that  was  observed  experimentally. 

4)  Ruppin  (1979)  calculates  the  absorption  due  to 
agglomerated  oxide  coated  A1  particles  using  the  Bruggeman 
theory  and  the  Maxwell-Gar  net t  theory.  For  small  particles 

O  # 

(50  A  in  diameter)  the  agreement  between  his  theory  for 
elongated  clusters  and  the  experimental  results  of  Granqvist 
et.al.  is  not  very  good.  For  larger  particles  (375  A  in 
diameter),  he  gets  good  agreement  with  the  experiment  when 
he  uses  the  Bruggeman  treatment  for  elongated  chainlike 
clusters  (with  the  ratio  of  the  long  axis  to  the  short  axis 
equal  to  30)  or  for  disclike  clusters. 

5)  In  a  later  paper,  Simanek  (1980)  proposes  that 
the  anomalously  high  absorption  is  due  to  Mott's  (1970)  a.c. 
conductivity  mechanism.  The  reason  for  this  new  model  is 
that  preliminary  measurements  by  Sievers  (1978)  indicated 
that  the  absorption  in  the  far  infrared  by  pure  amorphous 
aluminum  oxide  was  far  too  low.  In  this  new  calculation  the 
absorption  is  primarily  due  to  defects  (excess  aluminum  and 
oxygen  vacancies)  that  occur  at  the  interface  between  the 
aluminum  and  the  oxide.  Simanek  again  uses  the  model  of 
long  oxide  cylinders  with  spherical  metal  grains  mixed  in. 
With  the  use  of  Mott's  a.c.  conductivity  formulation  he 
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calculates  the  absorption  and  claims  to  get  good  agreement 
with  experiment. 

In  this  study  we  have  extended  the  measurements  of 
Tanner  et.al.  and  Granqvist  et.al.  by  measuring  the 
absorption  o£  small  aluminum  particles  from  the  far  infrared 
to  the  vacuum  ultraviolet  (.12  ym  to  250  ym) .  In  an  attempt 
to  resolve  the  discrepancy  between  the  theory  and  the 
experimental  results  use  is  made  of  a  simple  treatment  of 
nonspherical  shapes  (Huffman  and  Bohren  1980)  .  Although  the 
aluminum  particles  produced  are  known  to  be  almost  spherical 
(Granqvist  and  Buhrman  1976)  they  clump  together  to  form 
chains  and  clusters  which  behave  much  like  single  particles 
of  nonspherical  shape. 

Due  to  the  oxide  coating  on  the  particles  it  was 
necessary  to  consider  two  layer  particles  in  the  shape 
distribution  calculation.  This  requires  the  optical 
constants  of  the  amorphous  oxide  coating.  In  the  infrared 
these  were  determined  experimentally  by  doing  absorption 
measurements  and  the  results  of  Hagemann,  Gudat,  and  Kunz 
(1974)  were  used  in  the  ultraviolet. 


3 .  THEORY 


When  dealing  with  the  optical  properties  of  small 
particle  systems  it  is  important  to  understand  how  light 
interacts  with  bulk  solids.  This  chapter  begins  by 
introducing  the  terminology  used  when  dealing  with  optical 
properties  of  matter,  in  the  next  section  we  will  begin  the 
discussion  of  small  particles  by  looking  at  the  Mie  theory 
for  spheres  and  find  the  Rayleigh  limit  for  the  case  of  A1 
particles.  We  then  discuss  the  Rayleigh  theory  for 
ellipsoids  and  generalize  it  to  take  into  account  a 
distribution  of  ellipsoidal  shapes.  Finally,  this  shape 
distribution  is  used  to  calculate  the  absorption  of  small 
particles  of  Al203  and  comparison  is  made  with  experimental 
results . 

3.1  Optical  Properties  of  Bulk  Solids . 

Maxwell's  equations  for  the  interaction  of 
electromagnetic  radiation  with  a  nonmagnetic,  isotropic 
material  (in  cgs  units)  are: 
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(1) 


1  3  H 

C  'e  t 


£  5  £  +  4t rq 
c  5  t  c 


where 

E  is  the  electric  field  vector, 

H  is  the  magnetic  field  vector, 
c  is  the  speed  of  light  in  vacuo, 
e  is  the  dielectric  constant,  and 
o  is  the  electrical  conductivity. 

We  can  find  a  wave  equation  for  E  by  taking  the  curl  of 
and  substituting  (4)  into  the  resulting  expression, 


+  4ttct  5  E 


C.  3  t 


For  monochromatic  light,  of  wavector  q  and  frequency 
0) ,  a  single  plane  wave  can  be  used  to  describe  the  electric 


field  vector, 


E  *  Eq  exp  i(q*r-wt) 


Substituting  this  into  (5)  we  find, 


q2  ,  (e+  iiH, 


We  can  define  a  complex  index  of  refraction  m,  such  that, 

q  =  —  =  —  (n+ik)  A 

^  c  c 


where  n  is  just  the  usual  index  of  refraction,  k  is  normally 

called  the  extinction  coefficient,  and  q  is  a  unit  vector  in 

the  direction  of  q.  Substituting  this  in  the  wave  equation 

# 

(5),  the  expression  for  the  electric  field  becomes, 

E  «  Eq [exp- (  ~  k§* r  ) ]  [exp  i (  —  n§ • r-wt ) ]  ( 9 ) 

We  see  that  k  determines  the  extent  to  which  the  amplitude 
of  the  light  wave  is  attenuated  in  the  material,  i.e.  it  is 
a  measure  of  the  absorption  by  the  material.  c/n  is  the 
speed  of  light  in  the  material. 

Prom  (7) ,  we  see  that  a  complex  dielectric  constant 
can  be  defined, 

2  2 

q2  -  “T  e  -  ~r  (e'+ie")  (10) 

where'  e’  and  e"  are  the  real  and  imaginary  parts  of  c, 
respectively.  e'  is  just  the  old  t  from  Maxwell's 

equations,  and  e"s.*4it0/w.  By  squaring  (8)  and  comparing  to 
(10)  we  can  find  the  relationships  between  the  complex 
dielectric  constant  and  the  complex  index  of  refraction, 


e '  ■  n2-  k2 
e"  »  2nk 


(11) 


Either  of  these  two  sets  of  quantities  (e',c")  or  (n,k), 
characterize  the  optical  properties  of  the  material.  In 
discussions  of  optical  properties  of  small  particles  we 
usually  use  the  dielectric  constant,  e. 


3.2 


Introduction  to  Small  Particles  and  the  Mie  Theory. 


i. 


When  light  is  incident  on  a  system  of  small 
particles  it  can  either  be  scattered  or  absorbed.  The  sum 
is  called  extinction,  i.e., 

extinct ion* scat ter ing+absorpt ion . 

The  amount  of  light  transmitted  through  a  collection  of  n 
spherical  particles,  with  path  length  1  can  be  written, 

T  *  exp  (-NQextiTa2 )  1  (12) 

where 

N  is  the  number  of  particles  per  unit  volume*n/V  and 

a  is  the  radius  of  the  particle. 

Qext  is  called  the  efficiency  factor  for  extinction,  and  is 
equal  to  the  sum  of  the  absorption  efficiency  (Qabs)  and  the 
scattering  efficiency  (Qscat)*  These  quantities  will  be 
discussed  .later  in  more  detail.  It  is  more  convenient  to 

write  (12)  in  terms  of  the  mass  per  unit  area  of  particles, 

o,  where 

a  *  m/A  *  ml/V  *  (n4/3TTa 3  p)  1/V 


Here  p  is  the  bulk  density  of  the  material.  Therefor--1 ,  the 
transmission  equation  (12)  becomes, 


or 


T 


exp  (-3/4  ^ext  o) 
pa 


T  =  exp  (-aa/p) 


(13) 
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where  a=3Qext/4a  :''s  tbe  extinct^on  cross  section  per  unit 
volume  of  material. 

In  1908  Gustav  Mie  developed  a  theory  to  calculate 
^ext  an<^  Qscat  ^or  arbitrarily  sized  spherical  particles. 
The  derivation  is  given  many  places  (see  for  example  van  de 
Hulst  1957);  we  will  just  quote  the  results.  The 
extinction  cross  section  for  a  sphere  (where  Cext3  713  ®ext  ^ 


'ext  x‘ 


l  (2n+l)  Re  (a  +b  ) 
L  n  n 


and  the  scattering  cross  section  (cscat“7Ta2®scat^  *s 


'scat 


l  (2n+l)  (  |  a„  |  J  +  |  t>n  | 2 ) 


where  x  is  the  size  parameter  and  is  related  to  the 
wavelength  of  light  (A)  through, 

x  «  2ira/A  (16) 

The  Mie  coefficients  an  and  fc>n  are  complex  (and  complicated) 
expressions  involving  spherical  Bessel  and  Hankel  functions 
and  their  derivatives. 

It  turns  out  that  (14)  and  (15)  can  be  greatly 
simplified  if  the  following  two  conditions  hold, 
i )  x  <  <L , 
i  i )  |  m  |  x  <  <■  1 . 

The  first  condition  requires  that  the  size  of  the  particle 
be  small  compared  to  the  wavelength  of  light  and  the  second 
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condition  requires  the  phase  shift  of  light  within  the 

particle  to  be  negligible.  We  can  see  the  simplification  by 

looking  at  Fig.  2.  where  Cext/vol  (ct)  for  Al  vs.  particle 

radius  (a)  has  been  plotted  on  a  log-log  scale,  using  the 

Mie  theory,  for  several  different  wavelengths  of  light.  For 

each  of  the  wavelengths  there  is  a  flat  region,  where  the 

extinction  is  almost  independent  of  the  size  of  the 

0 

particle.  Such  a  flat  region  occurs  for  a<100A.  This 
region  is  where  the  above  two  conditions  are  fulfilled  and 
is  due  primarily  to  the  ai  term  in  the  Mie  theory.  It  is 
generally  called  the  Rayleigh  region,  after  Lord  Rayleigh, 
who  solved  the  problem  for  very  small  particles  many  years 
before  the  Mie  theory.  The  large  bump  for  the  three  curves 
with  X  >  lum  is  primarily  due  to  the  bt  term  in  the  Mie 
theory.  Physically,  it  is  caused  by  eddy  current  losses. 
Eddy  currents  are  important  when  the  diameter  of  the 
particle  is  close  to,  but  less  than,  the  skin  depth  of  the 
particle  material.  For  X  «,1  m  the  skin  depth  is  less  than 
the  radius  (for  a»100A)  so  that  the  interior  of  the  particle 
is  shielded  from  the  external  electric  field.  In  this  case, 
the  eddy  currents  are  confined  to  the  surface  so  that  the 
volume  of  material  able  to  absorb  energy  is  reduced. 
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Dependence  of  Absorption  on  Particle  Size 

Calculated  absorption  for  spherical  particles 
vs.  particle  radius  using  the  Mie  theory 
for  various  wavelengths  of  light 


It  can  be  shown  that  the  absorption  and  scattering 
cross  sections  calculated  in  the  Rayleigh  limit  are 
equivalent  to  considering  the  particle  as  an  electric  dipole 
in  a  uniform  electric  field  (see,  for  example,  van  de  Hulst 
1957) .  The  electric  polarizability  (a) ,  for  a  sphere 
embedded  in  a  medium  of  dielectric  constant  em,  is  (see 
Jackson  1962) , 
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where  e  is  the  complex  dielectric  constant  of  the  particle, 
and  a  is  the  radius  of  the  particle.  Cscat  and  cabs  are 
related  to  the  polarizability  through  the  relations  (van  de 


Hulst  1957)  , 
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where  k*  2 tt / X .  By  substituting  (17)  into  (18)  and  (19)  we 
get, 
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In  (20)  we  see  the  familiar  result  that  the  scattering  is 
proportional  to  l/X1*  (since  x=  2na/\),  provided  e  and  t 
are  slowly  varying. 
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3.3 


Rayleigh  Theory  for  Ellipsoids . 

So  far  in  our  discussion  of  small  particles  we  have 
been  dealing  with  spherical  particles.  Real  particles, 
however,  are  not  usually  spherical.  The  Mie  theory  is 
useful  for  calculations  for  spheres  but  is  quite  cumbersome 
and  not  easily  generalized  to  other  shapes  of  particles 
(except  infinite  cylinders)  .  As  discussed  in  the  previous 
section,  an  approximate  expression  for  the  absorption  and 
scattering  cross  sections  can  be  found  by  solving  the 
electrostatic  boundary  value  problem  to  obtain  the 
polarizability. 

We  can  generalize  this  to  other  shapes  of  particles, 
specifically  to  ellipsoids,  by  solving  the  boundary  value 
problem  in  ellipsoidal  coordinates.  This  derivation  is 
given  in  Stratton  (1941) .  The  expression  for  the 
polarizability  of  an  ellipsoid  with  semi-axes  a,b,c,  with 
the  field  parallel  to  the  j th  axis,  is, 
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where 


e  +L  ,  (e-e  ) 
m  .  j  m 


j«l,2,3  and 

V  is  the  volume  of  the  ellipsoid,  V=  4/3irabc. 

The  Lj's  are  called  depolarization  factors  and  are  related 
to  the  semi-axes  of  the  ellipsoid.  L3  ,  for  example,  is 
given  by, 
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They  are  normalized  such  that,  Lj+L2 +1,3*1.  Using  (18)  and 
(19)  we  can  get  expressions  for  C  and  Ca!;)g  for  a  single 
ellipsoid  with  arbitrary  Lj  , 
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where  the  superscript  j  means  the  light  polarization  is 
parallel  to  the  jth  axis. 

For  most  of  the  spectrum,  when  the  size  of  the 
particle  is  small  enough  the  extinction  is  almost  entirely 
due  to  absorption.  We  can  see  this  by  taking  the  ratio  of 


Cscat  to  Cabs ' 


Cabs  *  6ir 


( e ' -c " ) +e " 2 


where  the  real  and  imaginary  parts  of  c  are  explicitly 

shown.  Values  of  this  ratio  are  given  in  Table  1  using  the 

dielectric  constants  of  Al,  with  e  =1.0. 

Since  (25)  does  not  depend  on  the  shape  of  the 

particle  (i.e.,  there  is  no  Lj  dependence),  we  have  assumed 

0 

a  spherical  particle  with  a  radius  of  50A.  From  the  table 
we  see  that  it  isn't  until  well  into  the  UV  that  the 
scattering  becomes  important  and,  therefore,  in  general., 
scattering  may  be  neglected. 


Table  1 . 


Ratio  of  Scattering  to  Absorption 

for  Al  in  the  Rayleigh  Limit 


Wavelength  (vim) 

Frequency  (cm-1) 

Cscat//("abs 

100 

102 

10-7 

10 

103 

10-5 

1 

lO" 

10“  3 

.33 

3*10k 

10-2 

.10 

10  s 

.  3 

For  a  collection  of  randomly  oriented  identical 
ellipsoids  we  can  find  the  average  absorption  cross  section 
by  summing  (24)  over  j  and  dividing  by  three, 
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which  is  now  written  in  real  notation.  It  should  be  noted 
that  we  are  assuming  single  particle  absorption,  i.e.  the 
particle  separation  is  sufficient  that  the  field  seen  by  one 
particle  is  unaffected  by  the  presence  of  neighboring 
particles.  We  expect  the  absorption  to  be  large  whenever 
the  denominator  of  this  expression  is  close  to  zero.  There 
will  be  a  resonance  whenever  the  term  in  square  brackets  is 
equal  to  zero,  the  strength  of  the  resonance  depending  on 
the  size  of  e".  Setting  the  term  in  square  brackets  equal 
to  zero  and  solving  for  s' ,  we  get, 
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as  the  condition  for  resonance. 

For  a  sphere  all  of  the  Lj  '  s  are  equal,  and 
therefore  have  the  value  1/3,  so  there  will  be  a  maximum  in 
the  absorption  whenever, 

e '  »  -2e  . 

m 

For  a  distribution  of  identical,  noninteracting,  ellipsoids 
with  arbitrary  Lj ' s  there  will  be  three  peaks  in  the 
absorption.-  corresponding  to  the  three  values  of  Lj  . 
Furthermore,  we  see  that  for  a  system  of  particles  which 
contains  all  shapes  of  ellipsoids  (so  that  Lj  can  take  on 
any  value  from  0  to  1)  there  should  be  absorption  over  the 
entire  region,  where  e  '  is  negative. 

To  illustrate  these  three  cases  we  have  plotted  the 
absor ption/vol  for  Al  particles  in  air,  so  that  e  «1.0. 
(see  Fig.  3.  )  .  We  have  also  plotted  c  '  vs.  frequency  on 
the  same  frequency  scale  below.  For  the  spheres  curve  there 
is  a  large  peak  in  the  absorption  at  the  frequency  for  which 
e'=~2.  Physically,  this  corresponds  to  collective 
oscillations  of  the  free  electron  plasma  within  the 
boundaries  of  the  spherical  particles.  The  dashed  curve  is 
for  ellipsoids  with  L.=.01,  L2=.3,  and  L3=.69.  As  expected, 
there  are  three  pee  s.  A  discussion  of  how  the  remaining 
curve  was  obtained  is  left  for  the  next  section. 
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Calculated  Absorption  for  Al  in  the 
Rayleigh  Approximation 


Top  figure:  Calculated  Absorption  for 
spheres  (dot  dash),  a  single  ellipsoid 
with  Li=.01,  Lz=.3,  and  Lj=.59;  and  for 
a  continuous  distribution  of  ellipsoids 
(solid  line) 

Bottom  figure:  t'  vs.  frequency;  the  cross 
marks  the  point  where  c'=-2 


3.3.1  The  Shape  Distribution. 

The  average  absorption  cross  section  for  a 
collection  of  randomly  oriented,  identical  ellipsoids  is 
given  by  (26) .  Let  us  assume  that  the  small  particle  system 
under  consideration  consists  of  a  variety  of  shapes  of 
ellipsoids.  How  do  we  handle  this  situation?  If  a 
mathematical  expression  for  the  distribution  of  shapes  is 
known,  then  (26)  could  be  integrated  over  that  shape 
distribution  function,  i.e., 


«cab5  »  '  If  P<W  <Cab, 
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where  P  is  the  shape  distribution  function.  The  <<  >> 
denotes  the  double  averaging  that  is  being  done,  i.e.,  over 
orientation  and  shapes.  Because  of  the  restriction  that 
L  j  +L  2  +L  3  al ,  we  can  choose  L3  to  be  fixed  by  the  condition 
L3  =  1-Li-L2,  so  that  P  is  a  function  only  of  L  1  and  Lj. 

The  evaluation  of  this  integral  is  given  by  Huffman 
and  Bohren  (1980)  .  They  assumed  a  uniform  shape 
distribution  function  (i.e.,  P(LlfL2)*a  constant).  In  other 
words,  all  shapes  of  ellipsoids  are  taken  as  equally 
probable.  The  result  of  this  integration  is, 
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log (e/g  ) 
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In  the  log  term  the  principal  value  must  be  used,  i.e., 

!  C  I  >  0  ,  -  TT  <6-  <  TT  . 
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Let's  return  to  Fig. 


3. 


The  curve  labeled  CDE 


(which  stands  for  continuous  distribution  of  ellipsoids)  was 
.found  by  evaluating  (29)  .  There  is  absorption  ranging  from 
zero  frequency  (v»0)  ,  where  e'=  to  v*  1.2*105  cm-1  , 
where  c'  crosses  the  zero  axis.  The  maximum  absorption  for 
the  CDE  curve  is  approximately  an  order  of  magnitude  smaller 
than  the  maximum  for  spheres.  The  CDE  curve  has  a  broader 
absorption  band  though,  so  that  away  from  the  peak  it  has  a 
larger  absorption.  To  make  this  point  clearer  we  have 
replotted  this  graph  on  a  log-log  scale  (Fig.  4.  ).  At 
10um  (1000  cm-1)/  for  example,  the  absorption  calculated 
using  the  shape  distribution  is  approximately  three  orders 
of  magnitude  higher  than  for  spheres.  The  significance  of 
this  observation  will  become  clearer  in  Chapter  5. 

3.3.2  An  Example  using  the  Shape  Distribution. 


In  this  section  we  will  calculate  the  absorption  for 
oiAIjjOj  (corundum  crystal  structure)  particles  using  the 
continuous  shape  distribution  and  compare  it  to  experimental 
results.  The  aAl203  used  was  obtained  from  the  Buehler  Co. 
(alpha  micropolish) ,  the  particle  size  given  by  them  is 
.3um.  Measurements  were  done  using  the  standard  KB r  pellet 
technique  (see  Chapter  4)  in  the  region  of  the  infrared 
absorption  band.  The  result  of  these  measurements  is  shown 
in  Fig.  5.  For  comparison  we  have  also  plotted  the 
calculated  absorpt ion/voL  for  spheres  and  for  a  distribution 
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Fi9*  5.  Extinction  of  AI2O} 

Comparison  of  measured  extinction  (solid 
line)  of  Al203in  KBr  (e  *2.3)  with 
calculations  for  spheresm(Note :  a  =2.1* 
105cm"‘)  and  a  distribution  of  ellipsoidal 
shapes  (CDE)  in  the  Rayleigh  Approximation 
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of  ellipsoids  using  the  optical  constants  of  Barker  (1963) . 
If  we  compare  the  experimental  results  to  the  absorpt ion/vol 
calculated  for  spheres,  the  agreement  is  very  poor.  The 
locations,  strengths,  and  widths  of  the  peaks  differ 
substantially.  There  is  much  better  agreement  between  the 
experimental  results  and  the  CDE  curve.  The  widths  of  the 
calculated  and  experimental  bands  are  very  close  and  the 
measured  peak  magnitude  is  only  about  50%  higher  than  the 
CDE  curve.  The  experimental  curve  also  exhibits  more 
structure  than  the  CDE  curve.  Possible  reasons  for  the 
differences  in  these  two  curves  are  uncertainties  in  the 
optical  constants  used  to  find  the  CDE  curve  and  the  fact 
that  the  actual  particle  shapes  may  not  correspond  exactly 
to  a  uniform  continuous  distribution  of  ellipsoids.  The 
good  agreement  between  the  CDE  results  and  the  experimental 
results  gives  us  reason  to  believe  that  small  particle 
systems  of  other  materials  may  best  be  described  by  using  a 
shape  distribution. 


4,  EXPERIMENT 


To  measure  the  extinction  by  small  aluminum 
particles  two  slightly  different  sample  preparation  methods 
and  four  measuring  instruments  are  used.  In  this  chapter  we 
will  discuss  how  the  Al  and  Al203  samples  are  made  and 
measured.  First  we  will  describe  the  small  particle 
production,  followed  by  a  description  of  how  the  particles 
are  suspended  so  that  the  transmission  through  them  can  be 
measured.  In  order  to  determine  the  extinction  from 
transmission  measurements  it  is  necessary  to  find  a,  the 
mass  per  unit  area  of  particles.  We  will  describe  how  o  is 
determined.  Then  the  instruments  which  are  used  to  measure 
the  transmission  will  be  discussed.  Finally,  we  will 
discuss  the  final  step  in  the  sample  preparation  for  far 
infrared  measurements. 

4 . 1  Particle  Production . 

One  of  the  easiest  ways  to  make  very  small  particles 
is  by  using  the  technique  of  inert  gas  evaporation  (Kimoto 
et.al.  1963  and  Yatsuya,  Kasukabe,  and  Uyeda  1973).  The 
size  of  particles  that  is  produced  depends  on  several 
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factors:  the  atomic  weight  of  the  gas,  the  pressure  of  the 
gas,  the  evaporation  rate  of  the  material,  and  the  distance 
above  the  source  the  particle  are  collected  at.  In  general, 
the  smallest  particles  are  made  by  using  a  low  pressure  of 
helium  gas  at  a  slow  evaporation  rate  with  the  collection 
plate  as  close  as  possible  to  the  source. 

The  procedure  for  making  the  small  aluminum 
particles  is  as  follows.  First  the  chamber  is  evacuated  to 
2  x  10“5  torr  using  an  oil  diffusion  pump.  Purified 
aluminum,  obtained  from  the  Matheson,  Coleman,  and  Bell  Co., 
is  heated  in  a  tungsten  wire  filament  until  melted.  After 
the  filament  has  cooled,  the  system  is  purged  with  helium  to 
reduce  the  partial  pressure  of  oxygen,  and  then  reevacuated. 
Now  the  valve  between  the  pump  and  the  bell  jar  is  closed 
and  one  torr  of  helium  gas,  as  measured  by  a  thermocouple 
gauge,  is  introduced  into  the  chamber.  The  current  through 
the  filament  is  slowly  raised  until  the  smoke  of  small 
aluminum  particles  begins  to  be  produced,  A  metal  plate 
with  either  glass  slides,  polyethylene  or  LiF  substrates 
attached  to  it  is  positioned  10cm  above  the  filament  to 
collect  the  particles.  When  the  desired  density  of 
particles  is  achieved,  the  current  through  the  filament  is 
turned  off.  After  the  filament  is  allowed  to  cool,  air  is 
let  into  the  system  and  the  substrates  are  removed  from  the 
metal  plate. 
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The  particles  were  produced  under  the  same 
conditions  for  every  sample  so  that  ideally  the  particle 
size  should  be  the  same  for  all  samples.  Extensive  studies 
on  the  formation  of  A1  particles  produced  by  inert  gas 
evaporation  in  He  have  been  done  by  Yatsuya,  Kasukabe,  and 
Uyeda  (1973) .  They  determined  that  the  particles  are  almost 
spherical  and  measured  how  the  particle  size  depends  on  the 
pressure  of  the  gas  (for  a  given  filament  temperature  and 
distance  from  the  source  to  collection  plate)  .  By 
extrapolating  their  results  to  the  pressure  used  in  our 
experiment  we  determined  the  particle  diameter  to  be 
approximately  50A.  However,  the  spherical  particles  clump 
together  to  form  chains  and  clusters  which  changes  their 
effective  size  (and  shape) .  For  this  reason  we  do  not 
believe  it  would  be  useful  to  do  accurate  measurements  of 
particle  sizes  on  an  electron  microscope. 

Amorphous  aluminum  oxide  (aAljOj)  particles  can  be 
made  using  a  similar  technique.  In  this  case,  10  torr  of 
air  and  50  torr  of  He  were  admitted  into  the  system,  and  the 
smoke  was  made  in  the  same  manner  as  for  Al.  A  tantalum 
boat  was  used  rather  than  a  tungsten  filament  since  tungsten 
oxidizes  to  produce  W03  particles.  The  aALzOj  smoke  was 
collected  on  a  metal  plate  and  then  scraped  off  and  stored 
in  a  small  vial.  An  xray  powder  pattern  was  taken  on  the 
particles.  There  were  two  faint  and  indistinct  lines 
located  at  the  positions  where  the  two  brightest  lines  for 
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crystalline  yAl20j  occur  suggesting  that  the  material  is 
mainly  amorphous  with  perhaps  some  yAl203  mixed  in. 

We  also  did  measurements  using  commercially 
available  AlsOj  powder  obtained  from  the  Linde  Co..  The 
diameter  of  the  particles  is  given  as  .05pm.  This  material 
appears  to  be  primarily  yAljOj  based  on  the  analysis  of  the 
powder  xray  pattern  and  also  since  it  has  sharper  structure 
in  the  IR  than  the  smoke. 

4.2  Suspension  of  Particles . 

In  order  to  measure  the  amount  of  light  that  is 
absorbed  by  a  collection  of  particles  they  must  be  suspended 
in  some  manner.  The  easiest  way  is  to  collect  them  on  a 
transparent  substrate.  For  measurements  in  the  ultraviolet 
this  method  works  well  since  the  particles  can  be  kept 
reasonably  well  isolated  and  there  is  still  a  measurable 
amount  of  absorption.  We  use  this  method  for  A1  particles 
in  the  UV  ( . 12pm- . 45 pm)  with  LiF  as  a  substrate.  Because 

the  absorption  is  found  to  decrease  with  increasing 

wavelength  for  Al  particles,  the  density  of  particles  has  to 
be  increased.  Since  we  wish  to  keep  the  particles  isolated, 

the  best  way  to  do  this  is  to  stack  layers  of  particles  on 

top  of  each  other  with  a  transparent  substrate  between  each 
layer.  A  convenient  material  to  use  for  measurements  in  the 
visible  to  the  far  IR  is  polyethylene.  The  number  of  layers 
of  polyethylene  used  depends  on  how  strong  the  particle 
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absorption  is.  To  keep  the  layers  from  falling  apart,  and 
to  ensure  that  the  A1  particles  are  totally  surrounded  by 
polyethylene,  they  are  fused  together  by  placing  them 
between  microscope  slides  and  applying  gentle  hand  pressure 
while  heating  them  or»  a  hot  plate. 

For  the  A1203  smoke  and  powder  we  use  the  standard 
KBr  pellet  technique  (see,  £ or  example,  Miller  1972)  for  * 

measurements  from  the  visible  to  the  mid  IR.  A  small  amount 
(.2mg-.3mg)  of  the  Al203  is  mixed  with  . 5gm  of  IR  quality 
KBr  in  a  small  vial  which  has  a  steel  ball  bearing  in  it  and 
is  shaken  for  5  minutes  in  a  Crescent  "Wig-L-Bug"  dental 
amalgamator.  The  mixture  is  then  pressed  into  a  13mm 
transparent  pellet  by  using  a  special  die  and  applying  12 
tons  of  force.  For  measurements  on  Al20,  in  the  far  IR  the 
polyethylene  powder;  method  is  used.  The  material  is  mixed 
with  . lgm  of  polythethylene  and  poured  into  the  KBr  die. 

The  die  is  placed  on  a  hot  plate  until  the  powder  melts.  It 
is  allowed  to  cool,  and  then  the  polyethylene  disc  is  j 

carefully  removed  with  the  aid  of  a  razor  blade.  j 

4.3  Mass  Calibration .  '  J 

The  determination  of  the  mass  per  unit  area  (a)  for  } 

the  Al  particles  is  a  difficult  problem.  Since  o  is  very 
low  ( . 5ug/cmJ-5 gg/cm2 )  it  is  not  easy  to  measure  directly 

k 

using  an  analytical  balance.  To  solve  this  problem  we  used  \ 

I  the  following  method  to  determine  the  mass  of  a  sample 
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optically.  Al  smoke  was  collected  on  20  l"xl"  pieces  of 
glass  microscope  slides  (generally  with  higher  densities 
than  will  be  used  in  the  experiment)  .  The  transmission 
through  each  smoked  slide  was  measured  in  the  visible  using 
a  Cary  14  spectrophotometer.  The  slides  were  weighed  with 
the  smoke  on  and  reweighed  after  the  smoke  had  been  removed, 
the  difference  is  the  mass  of  the  Al.  Then  the  extinction 
cross  section  per  unit  volume  is  found,  at  some  wavelength 
(we  chose  5000A) ,  by  using  the  equation, 

T  =  exp(-ao/p) 

for  each  of  the  20  slides.  Since  the  absorption  depends 
somewhat  on  the  mass  per  unit  area, a  vs.  a  is  plotted  and  a 
least  squares  fit  is  done.  The  slope  was  found  to  be  10.6 
cm“3/gm  with  a  y-intercept  of  2.8*10*'  cm”  \  Therefore,  we 
find  for  the  typical  masses  we  will  use  in  the  experiment 
that  a  (5000A) =2 . 85xio ’em* 1 .  Now  that  a  is  known  at  this 
wavelength,  o  can  be  determined  for  an  arbitrary  sample  by 
measuring  the  transmission  through  it  and  inverting  (13). 
This  is  the  method  used  for  the  Al  particles  in  both  the 
polyethylene  sheet  method  and  the  LiF  method. 

For  the  AI2O3  samples  the  smoke  or  powder  was 
weighed  directly  or,  an  analytical  microbalance . 
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4.4  Measuring  Instruments . 


Four  different  measuring  instruments  were  required 
to  measure  the  absorption  by  Al  particles  from  . 12gm-250um. 
The  instrument  used  in  each  spectral  region  along  with  the 
type  of  sample  used  is  summarized  in  Table  2. 

In  the  far  infrared  (25um-25Qpm)  measurements  were 
done  on  a  Beckman  IR  11  spectrophotometer.  It  was  modified 
to  use  a  liquid  He  cooled  doped  germanium  bolometer  as  a 
detector  with  the  output  sent  to  a  PAR  model  28  lock-in 
amplifier.  The  output  from  the  lock-in  was  sent  to  a 
digital  voltmeter  or  a  chart  recorder.  Because  there  are 
reflectance  losses  at  the  polyethylene  surfaces  a  reference 
blank  was  made  in  the  same  manner  as  the  samples  using  clean 
polyethylene.  To  determine  the  absorption  by  the  Al 
particles,  the  ratio  of  the  transmission  through  the  sample 
to  the  transmission  through  the  blank  is  found. 

A  Perkin  Elmer  398  spectrophotometer  is  used  in  the 
mid  IR  region.  It  is  a  dual  beam  instrument;  the  sample  is 
placed  in  the  front  beam  and  the  reference  blank  is  placed 
in  the  back  beam.  An  optical  attenuator,  or  comb,  in  the 
reference  path  moves  to  keep  the  intensities  of  the  two 
beams  equal.  The  position  of  the  comb  indicates  the 
transmittance  of  the  sample  and  is  sent  electrically  to  the 
recorder.  This  is  useful  since  polyethylene  exhibit:; 
several  strong  sharp  absorption  bands  in  this  region  of  the 
spectrum  and  the  machine  can  usually  compensate  for  these. 
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For  the  near  IR  to  the  visible  region  we  used  a  Cary 
14  spectrophotometer.  It  is  also  a  dual  beam  instrument 
which  automatically  compensates  for  variations  in  the 
transmission  of  the  reference  sample. 

A  McPherson  model  235  Seya  Namioka  monochromater , 

with  a  Bausch  and  Lomb  grating  (radius=498 . 1mm, 
0 

l  ,  „."1500A,  600  lines/mm),  was  used  for  measurements  in 
the  UV,  with  a  hydrogen  discharge  lamp  as  a  light  source. 
The  lamp  produces  a  continuum  from  5000A  to  1600A  and  a 
multiline  spectrum  below  160oA.  The  detector  is  a 
photomultiplier  tube  that  had  been  coated  with  the 
fluorescent  material  sodium  salicylate,  with  the  output 
signal  going  to  a  picoammeter.  To  compensate  for  nonuniform 
light  levels  produced  by  the  hydrogen  lamp,  part  of  the 
light  that  exits  from  the  monochromater  is  deflected  to  a 
second  identical  phototube.  The  voltage  to  both  tubes  is 
supplied  by  a  single  programmable  high  voltage  power  supply. 
A  servomechanism  keeps  the  current  through  the  second 
phototube  constant  by  controlling  the  voltage  of  the  power 
supply,  thus  normalizing  the  output  of  the  first  phototube. 

The  experimental  apparatus  was  interfaced  to  an 
Imsai  8080  microcomputer  by  Ballart  (1980)  to  facilitate 
data  collection.  The  output  from  each  phototube  is  sent  to 
a  separate  picoammeter.  Tne  signals  from  the  two 
picoammeters  are  sent  through  a  multiplexer  to  an  analog  to 
digital  converter  and  then  to  a  parallel  port  of  the 
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microcomputer .  When  the  wavelength  is  being  scanned  a  pulse 
is  sent  from  the  McPherson  to  the  computer  every  5A  so  that 
the  computer  can  keep  track  of  the  wavelength.  To  find  the 
absorption  by  the  aluminum  particles  alone,  first  the 
transmission  vs.  wavelength  of  a  clean  part  of  the  LiP  is 
measured,  then  we  measure  the  transmission  through  the 
smoked  part,  and  finally  a  zero  level  scan  is  done  with  the 
first  phototube  pushed  aside  so  that  the  light  does  not  hit 
it.  The  zero  level  scan  is  necessary  because  the  signal 
from  the  picoammeter  depends  on  the  voltage  across  the 
phototube,  which  causes  the  zero  level  to  change  as  the 
wavelength  is  scanned.  The  computer  stores  the  data  from 
each  of  these  three  runs  and  then  performs  the  calculations 
to  find  the  absorption.  The  results  are  then  plotted  on  an 
X-Y  recorder  which  was  interfaced  to  the  computer  via  two 
digital  to  analog  converters.  A  block  diagram  of  the 
McPherson  electronics  is  given  in  Fig.  6. 

LiF  was  used  as  a  substrate  for  measurements  in  the 
UV  since  it  is  transparent  further  into  the  UV  than  any 
other  readily  available  material.  Measurements  were  not 
possible  below  1200&  because  the  LiF  became  too  absorbing. 
4.5  Far  IR  Samples . 

Thirty-six  samples  were  produced  by  the  polyethylene 
sheet  method.  The  mass  density  per  layer  of  the  samples 
ranged  from  .  Sug/cnr’/layer  to  4ug/cm2/iaYer  •  Transmission 
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measurements  were  done  on  each  sample  from  the  visible  to 
the  mid  IR  ( . 4 pm- 2 5 urn) .  We  used  (13)  to  calculate  the 
absorption  and  plotted  absorption  vs  frequency  on  a  log-log 
scale.  The  curves  were  then  compared  to  see  which  ones 
looked  most  alike.  The  reason  for  this  is  that  there  is 
some  difference  in  the  absorption  depending  on  how  densely 
the  particles  are  packed  onto  the  polyethylene.  They  broke 
up  into  four  groups  (which  we  call  I-IV) ,  the  samples  with 
the  lowest  mass  densities  ( . 5ug/cm2/layer  to 
1 . 2yg/cm2/layer )  generally  fell  into  group  I  and  those  with 
the  highest  generally  fell  into  group  IV  ( 2yg/cm2/layer  to 
4yg/cm2/layer ) .  The  samples  in  group  I  were  then  fused 
together  using  heat  and  pressure  to  form  sample  I.  The  same 
thing  was  also  done  with  groups  II-IV.  These  samples  now 
have  a  hiyh  enough  mass  density  to  produce  appreciable 
absorption  in  the  far  infrared. 
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5.  RESULTS  AND  DISCUSSION 


5.1  A1  Results . 

In  Fig.  7,  we  have  plotted  the  result  of  the 
extinction  measurements  for  sample  I  for  the  wavelength 
region  .4pm-170pm  and  have  compared  it  to  the  sphere  theory 
and  to  the  shape  distribution  calculation  (CDE) .  This 
sample  has  an  average  mass/area/layer  of  . 7pg/cm 2/layer . 
Also  plotted  on  this  figure  is  the  average  extinction 
coefficient  for  five  samples  with  mass  densities  ranging 
from  lpg/cm2  to  5pg/cm2  in  the  wavelength  region  of 
.12pm~.5pm.  The  statistical  uncertainty  in  the  measurements 
is  smaller  than  the  size  of  the  data  points  for  50um<A < . 12pm 
and  increases  in  the  far  infrared  as  shown  by  the  error 
bars.  The  large  error  in  the  far  IR  is  caused  by  a 
relatively  large  electronic  noise  level.  Also  in  this 
region  the  A1  is  becoming  almost  transparent  so  the 

transmission  through  the  sample  is  very  close  to  the 

transmission  through  the  reference.  The  discontinuity  in 

the  data  at  .5pm  is  due  to  the  fact  that  for  shorter 
wavelengths  the  dielectric  constant  of  the  medium  is  1 

rather  than  2.3  (i.e.  we  switched  from  the  polyethylene 
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Fig.  7.  Measured  Extinction  for  Al  (sample  I) 

Comparison  of  experimental  results 
(open  circles, e  «2.3;  x's,e  -1)  to 
sphere  theory  (dashed  line)  and  CDE 
theory  (solid  and  dotted  lines) 


samples  to  the  LiF  substrate  samples.)  The  agreement  between 
the  data  and  the  CDE  curve  is  quite  good  and  is  certainly 
much  better  than  between  the  data  and  the  sphere  theory, 
particularly  in  the  far  IR  where  the  sphere  theory  is 
approximately  three  orders  of  magnitude  below  the  other  two 
curves.  The  most  obvious  difference  between  the  data  and 
the  CDE  curve  is  the  bump  at  12pm.  This  is  almost  certainly 
due  to  an  amorphous  Al20j  layer  on  the  particles  since 
aAl20j  has  an  absorption  band  in  this  region.  It  is  well 
known  that  A1  surfaces  oxidize  rapidly  and  several 
monolayers  of  oxide  can  form  on  a  surface  in  a  short  period 
of  time  even  in  a  fairly  good  vacuum. 

Fig.  8.  shows  the  measured  extinction  coefficient 
for  sample  IV,  which  has  c  »2 . 3pg/cm2/layer ,  in  the 
wavelength  range  .  33um-280uin.  In  general  sample  IV  has  a 
higher  extinction  coefficient  than  sample  I.  The  reason  for 
this  may  be  that  the  density  of  particles  is  so  high  that 
the  independent  particle  assumption  is  being  violated,  i.e., 
when  a  particle  goes  through  resonance  it  affects  its 
neighbors.  Another  possible  explanation  is  that  since  the 
particles  are  so  packed  together  they  effectively  form 
particles  of  a  larger  size  so  the  use  of  Rayleigh  theory  is 
no  longer  valid. 

In  future  discussions  in  this  work  we  will  only 
consider  sample  I  since  it  has  the  lowest  particle  density 
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and,  therefore,  will  be  most  likely  to  fulfill  the  Rayleigh 
conditions  and  the  noninteracting  particle  assumption. 

5.2  a12°3  Results* 

The  measured  extinction  cross  section/volume  for 
aAl20j  is  plotted  in  Fig.  9.  For  comparison  we  have  also 
plotted  the  measured  extinction  coefficient  for  ya120j.  The 
two  curves  are  almost  identical  from  .  6um  to  25ym  and  begin 
to  differ  in  the  far  IR,  with  the  amorphous  material  having 
a  higher  extinction  than  the  yAl20j.  This  is  consistent 
with  the  results  of  Strom  et.al.  (1974)  .  They  observed 
that  amorphous  materials  tend  to  have  larger  extinction  in 
the  far  IR  than  their  crystalline  counterparts;  even  orders 
of  magnitude  larger.  The  extinction  spectrum  for  aAl203  is 
a  composite  of  measurements  for  several  different  samples. 
For  the  region  of  the  IR  absorption  band  of  aAl203 
(10um-25um) ,  measurements  were  done  on  three  samples  with 
masses  ranging  from  ,207mg-. 345mg  with  KBr  as  a  medium. 
Since  the  extinction  decreases  significantly  on  both  sides 
of  this  band  it  is  necessary  to  use  a  higher  mass  of 
particles  to  produce  a  measurable  amount  of  extinction.  For 
measurements  in  the  wavelength  regions  25um-160pm  and 
.6um-10um  we  used  two  samples,  both  with  4.4mg  of  aAljOj,  in 
a  polyethylene  medium. 

From  Fig.  9.  we  see  that  the  extinction  has  a 
minimum  at  approximately  4pm  and  increases  for  shorter 
wavelengths.  Calculations  of  the  absorption  coefficient  for 
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aAljOj  (for  particles  of  any  shape)  show  that  tne  absorption 


should  continue  to  decrease,  so  that  the  increase  is 
probably  attributable  to  scattering.  Prom  calculations  of 
single  particle  scattering  it  is  found  that  the  scattering 
coefficient  begins  to  increase  at  about  Sum,  but  is  several 
orders  of  magnitude  lower  than  the  measurements.  This 
difference  in  magnitude  may  be  due  to  multiple  scattering 
since  there  is  a  fairly  high  concentration  of  particles 
(4.2%  by  mass).  In  the  region  of  the  minimum  the  total 
extinction  is  smaller  than  the  uncertainty  in  the 
measurements,  as  shown  by  the  error  bars. 

5.3  Two  Layer  Particles . 

As  discussed  earlier,  the  major  difference  between 
the  experimental  results  and  the  CDE  curve  is  the  oxide  bump 
at  12um.  Since  it  is  not  possible  to  avoid  the  oxide 
coating  we  have  decided  to  take  this  into  account  in  the 
shape  distribution  calculation.  It  should  be  noted  that  the 
absorption  of  a  coated  particle  is  not  simply  the  sum  of  the 
absorptions  by  the  core  and  the  coating,  since  the  electric 
field  seen  by  the  core  is  affected  by  the  presence  of  the 
coating.  Therefore,  we  consider  the  problem  of  a  coated 
ellipsoid,  i.e.,  an  ellipsoid  of  material  1  surrounded  by  a 
layer  of  material  2,  The  absorption  cross  section/volume 
for  a  collection  of  randomly  oriented,  noninteracting, 
identical  ellipsoids  is  (Bohren  and  Huffman,  to  be 
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where 


k  is  the  wavenumber  {k»2Tt/X)  , 

ej  is  the  complex  dielectric  constant  of  the  inner 

material , 

e2  is  the  complex  dielectric  constant  of  the  outer 

material , 

c  is  the  dielectric  constant  of  the  medium,  and 
in 

f  is  the  volume  fraction  occupied  by  the  inner 
ellipsoid . 


To  find  the  absorption  coefficient  for  a  collection 
of  randomly  shaped  ellipsoids  we  follow  the  same  procedure 
used  in  Chap.  3  ,  i.e.,  we  integrate  (30)  over  a  shape 
distribution  function.  Again  it  is  assumed  that  the  shape 
distribution  function  is  a  constant,  and  the  integration  is 
done  over  Lt  and  L2  (with  Lj  fixed  by  the  condition 
L  j»l-L 1 ~L 2 ) .  Because  of  the  complexity  of  (30)  it  was 
necessary  to  do  the  integration  numerically.  The  domain  of 
integration  is  shown  in  Fig.  10a.  Although  all  shapes  of 
ellipsoids  are  contained  in  any  of  the  regions  1-6,  for  ease 
of  computation  we  integrate  over  the  whole  triangle  (with 
proper  normalization).  The  double  integral  over  Li  and 
is  replaced  by  the  double  sum, 


51 


Region  of  Integration 


a)  All  shapes  of  ellipsoids  are 
contained  in  each  of  the 

six  regions 

b)  The  shaded  area  is  the  domai 
of  integration 
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with  increments  of  1/N,  where  N  is  the  number  of  steps  that 
Lj  (or  I* 2  or  Lj)  is  divided  into. 

Optical  constants  for  the  amorphous  oxide  coating 
and  the  aluminum  core  are  required  for  the  computations. 

For  Al,  a  Drude  free  electron  model  is  fitted  to 

experimental  data  tabulated  by  Hagemann  et.al.  (1974) .  The 
Drude  expression  for  the  complex  dielectric  constant  is 
(see,  for  example,  Wooten  1972), 


1  wu+irb) 

where  o  and  ru  are  the  plasma  frequency  and  the  damping 

p  D 

constant,  respectively.  We  achieved  a  good  fit  to  the 

measured  reflectance  of  Al  over  a  broad  energy  range  from  0 

to  15  eV  by  choosing  w  *15eV  and  I\a.6eV. 

P 

For  the  aAl203  coating,  the  dielectric  constants 


tabulated  by  Hagemann  et.al 


(1974)  are  used  for 


wavelengths  less  than  ,3um.  Experimentally  determined 
dielectric  constants  for  aAl203  do  not  appear  to  be 
available  in  the  infrared.  Since  bulk  samples  of  aAl203  are 
not  available  for  optical  constants  determinations  by  the 
usual  techniques  of  reflection  and/or  transmission 
spectroscopy,  we  were  forced  to  use  the  small  particle 
samples  described  in  Chap.  4.  Details  of  how  the  optical 
constants  were  determined  are  given  in  Appendix  A. 
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We  use  these  optical  constants  for  aAl203  and  Al  in 
(30)  and  (31)  in  order  to  do  a  numerical  calculation  for  the 
absorption  by  a  shape  distribution  of  coated  ellipsoids, 
with  N  and  f  as  adjustable  parameters.  In  Fig.  n.  we  see 
the  effect  changing  N  has  on  the  calculated  absorption  for 
Al  (f*1.0).  In  general,  the  numerical  result  follows  the 
analytic  curve  (Fig.  ll.d  )  until  some  critical  frequency 
where  the  slope  changes,  so  that  the  absorption  is  quadratic 
in  the  frequency  rather  than  linear.  Also  plotted  in  this 
figure  are  the  results  of  Granqvist  et.al.  (1976).  It  is 
seen  that  their  measured  far  IR  extinction  has  quadratic 
frequency  dependence.  Initially,  N  is  chosen  large  enough 
so  that  the  cutoff  occurs  past  the  oxide  absorption  band.  f 
is  then  varied  until  the  best  fit  to  our  experimental 
results  in  the  region  of  the  oxide  peak  is  obtained.  In 
this  way  it  was  determined  that  f«.5,  i.e.,  half  of  the 
total  particle  volume  is  occupied  by  aAl203  ,  Finally,  N  is 
adjusted  until  the  numerical  result  closely  fits  our  far  IR 
data  for  oxide  coated  Al  particles,  and  is  found  to  be  aoout 
800. 

The  limitation  on  the  size  of  N  that  was  used  may 
have  a  physical  interpretation.  Since  very  small  values  of 
L  are  cut  out  (see  Fig.  10. b  )  this  effectively  eliminates 
extremely  elongated  ellipsoids  (as  well  as  very  large  flat 
ellipsoids) .  In  fact  we  can  get  an  idea  of  the  maximum 
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Effect  of  Changing  N  on  Calculated 
Absorption 

a)  N*50  (dotted  line) 

b)  N=200  (dashed  line) 

c)  N  =  8  0  0  (dot  dash  line) 

d)  analytic  result 

Also  shown  is  the  measured  extinction  of 
Granqvist  gt  gl .  (1976)  (open  circles) 


length  of  ellipsoid  that  is  possible.  As  a  special  case 

let's  consider  a  prolate  spheroid  with  b*c  a  where  a,b,  and 
c  are  the  semi-axes  of  the  ellipsoid.  From  van  de  Hulst 
(1957),  the  expression  for  the  depolarization  factor,  Li,  is 

L  -  [-1+  ln<y±|)];  e2  =  l-(b/a)2. 
i  e2  2e  1-e 

For  Li  we  will  use  the  minimum  value  possible  1/2400.  By 
trial  and  error  we  find  a-lOOb.  Therefore,  the  maximum 
length  of  ellipsoid  possible  in  our  model  is  one  in  which 
the  longest  axis  is  100  times  the  smallest  axis. 

Fig.  12.  shows  the  final  fit  of  the  distribution  of 
coated  ellipsoids  to  the  experimental  results.  This  is,  in 
a  sense,  a  comprehensive  summary  of  this  work  from  which  we 
can  draw  reasonable  insights  concerning  the  long  standing 
problem  of  anomalous  far  IR  absorption  by  Al  particles.  The 
agreement  between  the  data  and  the  numerical  fit  is  seen  to 
be  very  good.  The  largest  discrepancy  occurs  in  the  UV  with 
a  maximum  relative  error  of  about  a  factor  of  2  (when 
compared  to  the  CDE  curve  with  cm*l) ,  but  at  these  short 
wavelengths  we  expect  the  assumptions  of  the  Rayleigh 
approximation  to  be  violated  (Fig.  2.  )  leading  to  an 
expected  decrease  of  volume  normalized  extinction  compared 
to  the  small  particle  limit  (Huffman  1977).  In  addition  the 
calculation  ignores  scattering  which  becomes  important  in 
the  UV .  This  leads  to  an  increase  in  the  extinction.  The 
combination  of  these  two  effects  probably  explains  the 


absence  of  the  predicted  dip  in  the  UV.  In  the  middle  IR 
the  aqreement  is  quite  good.  The  far  infrared  anomaly  is 
resolved  without  the  need  for  anything  other  than  a 
classical  shape  distribution  calculation.  There  was  no 
evidence  of  periodic  structure  in  the  far  IR  caused  by 
quantum  size  effects.  Of  course  we  would  expect  any  such 
structure  to  be  masked  by  the  large  absorption  due  to  shape 
effects , 

Some  objections  can  be  raised  to  this  simple 
explanation,  however.  After  all,  we  have  used  a  rather 
simplified  model,  i.e.  real  particles  are  not  ideal 
ellipsoids.  However,  there  is  independent  evidence  that  the 
absorption  determined  by  integration  over  a  wide  range  of 
ellipsoidal  shapes  closely  approximates  the  absorption  by  a 
collection  of  irregularly  shaped  particles.  The 
measurements  on  A120j  given  in  Chap.  3  ,  as  well  as 
measurements  on  crystalline  quartz  (Huffman  and  Bohren  1980) 
agree  well  with  the  CDE  model.  Another  possible  shortcoming 
of  our  model  is  that  we  treat  all  possible  shapes  of 
ellipsoids  as  being  equally  likely.  This  works  well  for  Al, 
aAl2Oj,  and  Si02,  but  it  may  be  that  a  nonuniform  shape 
distribution  function  will  have  to  be  used  for  some  other 
mater i als . 

One  thing  that  has  not  been  included  in  our  model  is 

the  effect  on  the  absorption  due  to  a  size  limited  mean  free 

path,  since  the  electron  mean  free  path  is  larger  than  the 
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size  of  the  particle  (see  Kreibig  1974) .  This  can  be  taken 


into  account  by  replacing  the  bulk  damping  constant  rb  in 


(32)  by , 


rb+vF/(Lc) 


where  vp  is  the  Fermi  velocity,  c  is  the  speed  of  light,  and 
L  is  a  geometrical  factor  proportional  to  the  size  of  the 
particle  (L»4/3a  for  a  sphere  of  radius  a).  The  complex 
dielectric  constant  now  becomes, 


oj [ui+i (rb+vp/ (lc)  ) ) 

To  find  the  effect  of  the  size  limited  mean  free  path  on  the 
absorption  we  substitute  (34)  into  (31)  for  cj.  To  estimate 
the  magnitude  of  this  effect  we  assume  a  spherical  particle 
of  radius  25A.  The  Fermi  velocity  for  Al  is  2*10e  cm/sec  so 


that  ro«  5Pb.  With  this  value  of  Pb 


the  absorption 


coefficient  at  lOOum  is  about  a  factor  of  2  larger  and 
decreases  at  shorter  wavelengths.  Since  the  original 
discrepancy  is  over  3  orders  of  magnitude,  this  factor  can 
be  ignored. 

The  mean  free  path  limitation  introduces  a  size 
dependence  into  the  expression  for  the  absorption.  For 
spherical  particles  small  enough  to  satisfy  the  Rayleigh 
criteria,  Granqvist.  (1978)  calculates  that  the  size 
dependence  of  P  leads  to  the  absorption  being  inversely 
proportional  to  the  size.  In  contrast,  the  absorption 


V  .*  •»  !  .«  ,» 


coefficient  for  pure  A1  particles  calculated  using  the  shape 
distribution  is  very  insensitive  to  the  size  of  r.  in  fact 
changing  r  by  a  factor  of  5  produces  virtually  no  change  in 
the  absorption  coefficient.  It  may  be  concluded  that  the 
absorption  coefficient  for  a  distribution  of  shapes  is 
relatively  size  independent.  This,  of  course,  assumes  that 
the  particles  are  small  enough  that  Rayleigh  theory  holds  so 
that  eddy  current  losses  are  negligible. 

5.4  Other  Explanations  of  "Anomalous  Absorption. 

As  pointed  out  in  Chap.  2  there  have  been  a  number 
of  theories  proposed  to  account  for  the  anomalously  high  far 
IR  absorption  observed  by  Granqvist  et.al.  (1976).  We  will 
now  discuss  these  in  further  detail  to  see  how  well  they 
agree  with  the  extension  of  the  measurements  to  shorter 
wavelengths  that  we  have  done.  Since  the  theories  all 
assume  the  dielectric  constant  of  the  medium  to  be  1,  we 
will  compare  their  results  to  the  experimental  fit  curve 
with  em«l  (rather  than  2.3). 

1)  In  Simanek's  (1977)  first  theory  he  proposed  that 
the  absorption  was  taking  place  mainly  in  the  oxide  coating 
of  the  particles.  He  modeled  the  particle  system  as  a 
collection  of  long  randomly  oriented  aAlaOj  cylinders  with 
spherical  Al  particles  embedded  within.  For  the  absorption 
coefficient  he  finds, 
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(35) 


a  (id)  -  ^1+ - - - lira  e 

(Re  £MG+l)2+(Im  eMG) 2 


where  s\._  is  the  complex  dielectric  constant  of  the 
MG 


composite  using  the  Maxwell-Gar  nett  theory.  eMr  is  given 


e  (l+2f)+2  (1-f) 


where  f  is  the  fraction  of  the  total  volume  of  the  cylinder 


occupied  by  the  metallic  cores.  In  the  far  infrared 


Ci  | >>  |  ca| (ei  is  the  complex  dielectric  constant  of  Al  and  e2 


the  complex  dielectric  constant  of  aAl20j>  and  can  be 

MG 


simplified , 


.  l  +  2f 


■  E  „  “  c  - R 


Because  we  do  not  want  to  limit  the  calculation  to  the  far 


IR  we  do  not  make  this  assumption,  but  rather,  use  the  full 


Maxwell  Garnett  dielectric  constant.  For  e2  the  dielectric 


constants  that  were  determined  experimentally  were  used 


(Appendix  A)  instead  of  the  approximate  far  IR  dielectric 


constant  used  by  Simanek, 


c  -  10  +  12X10"1  \ 
2 


The  result  of  this  calculation  is  shown  in  Fig.  13. 


comparison  we  have  also  plotted  the  result  of  the  numerical 


fit  with  c  *1.  In  the  far  IR  the  agreement  is  good  but  for 


shorter  wavelengths  it  is  very  poor.  We  can  conclude  that 
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Fig.  13.  Simanek's  (1977)  Results 

Simanek's  results  (solid  line)  compared 
to  experimental  fit  with  e  «1 
(dashed  line)  n 
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the  oxide  coating  is  not  the  major  contributor  to  the 
absorption. 

2)  Glick  and  Yorke  (1977)  did  a  calculation  to 
determine  the  lattice  contribution  to  the  absorption  for 
small  Al  particles  since  they  believed  the  electronic 
component  to  be  small.  Their  results  showed  that  the 
lattice  absorption  should  mirror  the  phonon  density  of 
states  for  Al.  In  Fig.  14. a  the  density  of  states  for  Al 
is  plotted.  We  have  plotted  their  calculated  absorption  (as 
far  as  it  extends)  in  Fig.  14. b  and  have  compared  it  to  our 
experimental  fit  curve.  Their  result  is  about  a  factor  of  6 
higher  than  the  experimental  data  and  does  not  have  quite 
the  same  frequency  dependence.  They  suggest  that  further 
measurements  be  done  to  see  if  there  is  a  leveling  off  in 
the  absorption  that  they  predict  since  the  data  of  Granqvist 
et.al.  stops  just  short  of  this  point.  We  have  extended 
the  measurements  to  shorter  wavelengths  and  do  not  see  a 
leveling  off.  Furthermore,  it  would  seem  that  if  Glick  and 
Yorke  extended  their  calculations  the  absorption  should 
begin  to  decrease  significantly.  This  was  also  not  observed 
experimentally  (compare  Fiy.  7.  ).  Due  to  this  discrepancy 
it  appears  that  the  theory  of  Glick  and  Yorke  does  not 
adequately  represent  the  absorption  by  small  Al  particles. 

There  are  some  theoretical  reasons  to  question  GLick 
and  Yorke' s  result.  For  one  thing  they  used  the  bulk 
density  of  states  rather  than  the  actual  small  particle  mode 


LOG  FREQUENCY  (cm*') 

The  Theoretical  Results  of  G.lick 
and  Yorke 

a)  Density  of  states  for  M 

b)  Click  and  Yorke's  results  (solid 

line)  compared  to  our  experimental 
fit  with  -  *1  (dashed  line) 


density  in  their  calculation.  Simanek  (1980)  has  pointed 


out 

that  for  very 

small 

particles 

there  is  another 

term 

in 

the 

vibrational 

mode 

density 

due 

to  surface 

modes . 

He 

concludes  that  for 

particles  with 

a 

diameter  of 

50A 

the 

surface  term  is  five  times  as  yreat  as  the  volume  term 
considered  by  Glick  and  Yorke.  Also  this  surface  term  is 
linear  in  frequency  which  implies  a  strong  deviation  from 
the  quadratic  dependence  seen  experimentally. 

Another  apparent  problem  with  Glick  and  Yorke' s 
calculation  is  that  they  assumed  the  electric  field 
penetrates  only  a  few  angstroms  into  the  particle.  This 
assumption  is  based  on  a  paper  by  Rice,  Schneider,  and 
Strassler.  They  calculated  the  effect  of  reducing  the 
particle  size  on  the  polarizability  and  found  that  reducing 
the  size  decreases  the  polarizability,  thus  allowing  the 
electric  field  to  penetrate  the  particle.  Their 
calculation,  however,  assumes  a  perfect  conductor. 
Therefore,  it  is  questionable  whether  the  penetration  depth 
calculated  in  this  manner  can  be  properly  used  for  Al . 

3)  Lushnikov  et.  al.  (1978)  calculated  the 
absorption  by  small  particles  of  Al  by  taking  into  account 
the  Coulomb  interactions  between  conduction  electrons. 
Their  calculation  gives  quadratic  frequency  dependence  of 
the  absorption  but  is  two  orders  of  magnitude  smaller  than 
the  experimental  results.  They  admit  that  the  model  they 
use  can  not  be  expected  to  give  good  quantitative  agreement 


with  experimental  data.  They  say  that  this  is  because  their 
model  assumes  spherical  pure  A1  particles  of  a  definite 
size,  whereas  real  particles  have  an  oxide  coating  and  tend 
to  clump  together  into  clusters  and  chains. 

4)  Ruppin  (1979)  extends  the  theory  of  Simanek 
(1977)  to  various  metal-oxide  configurations  by  considering 
the  particles  to  cluster  in  the  shapes  of  prolate  spheroids. 
Two  models  are  used  for  the  topology.  In  the  first,  the 
metal  particles  are  embedded  in  an  amorphous  oxide  matrix 
(Maxwell-Garnett  theory).  The  second  model  consists  of 
metal  and  oxide  particles  forming  clusters  (Bruggeman 
theory).  To  calculate  the  absorption  he  uses  (26),  with 

“.bs  ■<Cabs>''v- 

aabs  *  5 

j-1 

where  e„  has  been  taken  to  be  1.  c  is  the  complex 

dielectric  constant  of  the  aggregate  which  for  the 
Maxwell-Garnett  theory  is  equal  to  and  is  given  by  (36). 

MG 

For  the  Bruggeman  theory  ea  ■  and  is  determined  by  the 
equation , 

3f  ,  3  (1-f ) 

■  1 

where,  as  usual,  ci  and  e2  are  the  complex  dielectric 

constants  of  the  metal  and  oxide  respectively.  Since  the 

dielectic  constants  of  aAl20j  were  not  known,  Ruppin  used 

the  approximate  dielectric  constant  from  Simanek  (1977), 


given  in  (38).  He  suggests,  however,  that  this  dielectric 
constant  is  a  possible  source  of  error  in  his  calculation. 
The  results  of  Ruppin's  calculations  are  given  in  Fig.  15. 
(with  the  ratio  of  the  major  to  the  minor  axis  equal  to  6) 
where  we  have  used  the  dielectric  constants  that  we 
determined  experimentally  for  aAl20i.  For  comparison  the 
experimental  fit  curve  to  our  data  with  e  «1.0  is  also 
given.  In  the  far  IR  there  is  reasonable  agreement  but  for 
shorter  wavelengths  both  of  his  curves  differ  significantly 
from  the  experimental  results,.  From  this  we  conclude  that 
the  absorption  calculated  assuming  a  collection  of  identical 
ellipsoids  moti"!i*d  by  either  the  Maxwell-Garnett  or  the 
Bruggeman  theories  does  not  adequately  explain  the 
experimental  results. 

5)  In  Simanek's  second  theory  (1980)  the  absorption 
is  calculated  using  Mott's  (1970)  a.c.  conductivity 
formulation.  He  abandons  his  previous  theory  because 
preliminary  measurements  by  Sievers  (1978)  suggest  that  the 
absorption  by  pure  a<\l2Oj  particles  is  far  too  low.  Our 
measurements  show,  however,  that  aAi203  particles  have  an 
absorption  coefficient  comparable  to  that  of  the  oxide 
coated  Al  particles  in  the  far  IR.  In  Simanek's  new  model 
the  mechanism  for  absorption  involves  electrons  which  are 
localized  at  defects  of  the  oxide  coating.  These  defects 
are  due  to  excess  aluminum  and  oxygen  vacancies  which  arise 
because  of  nonstoichiometry  of  the  aAl;0,  at  the  metal-oxide 
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interface.  The  absorption  is  caused  by  electrons  executing 
resonant  transfer  between  trap  states  localized  at  the 
defect  sites.  The  a.c.  conductivity  for  this  mechanism  has 
been  found  by  Mott  to  be, 

Re  cr(uj)  -  3TTe2/ft  [N  (EpHSw]  2  a-  5  [  In  }  4  (39) 

where 

1/a  is  the  range  of  the  electron  wave  function, 

I0  is  related  to  the  overlap  of  localized  electron  wave 
functions , 

N(Ep)  is  the  density  of  trap  states  at  the  Fermi  energy 
and  is  assumed  to  be  constant. 

The  absorption  coefficient  is  given  by  (35) ,  where 
for  eMG  we  have  used  (37)  which  is  valid  for  wavelengths 
less  than  Sum. 

a(w)  =  t  R  Im  e  (1+  - - - ]  .  (40) 

J  2  (R  Ro  e  +i) 2+ (Im  e  ) 2R2 

?  j 

Simanek  used  (35)  to  calculate  the  far  IR  absorption  with 
lrac2  replaced  by  4ffa/<jj  ,  and  the  following  assumptions, 

1)  1/a  “  2A, 

2)  R*1.7,  which  corresponds  to  f=.2, 

3)  (R  Im  C2)2  in  the  denominator  of  the  bracketed 
expression  is  small  and  Re  a  is  taken  to  be  10, 

4)  I  o  =  10  eV. 

He  then  determined  N(e  )  to  be  2 <1 0 2  1  cm" 5eV~ 1  by  requiring 

r 
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agreement  with  the  data  ok  Granqvist  et.al.  (1976).  He 
states  that  this  value  is  typical  for  metal-oxide 
interfaces,  citing  the  work  of  Halfcu'itter  (1974)  wno 
considered  Nb-NbiOs  interfaces.  We  have  found,  however t 
that  by  using  Simanek's  values  l)-4)  that  N(Ep)  should  be 
3*10* °cm“ 3eV“ 1  to  fit  Granqvist's  data  at  lOOym. 

To  extend  Simanek's  (1980)  theory  to  shorter 
wavelengths  the  Ime2raust  include  a  term  to  account  for  the 
dielectric  constant  of  the  pure  aAl2Oj , c2pure  ,  i.e., 

Im  e  “  Im  c  ^5Ure-*-4‘rra/uj , 

2  2 

With  this  expression  for  Im  in  (40)  we  calculated  the 
absorption,  which  is  plotted  in  Fig.  16  .  We  took  N(Ep) 
to  be  3*102 °cm" JeV" 1  and  used  our  experimentally 
determined  dielectric  constants  for  Re  e?  and  lmc2PUre,  along 
with  Simanek's  values  for  1/a  and  l0  and  R.  For  comparison 
our  experimental  fit  curve  with  e  »1  is  also  shown.  We  see 
that  there  is  reasonably  good  agreement,  the  primary 
difference  being  the  location  of  the  oxide  bump.  A  more 
serious  problem  with  this  model  is  its  strong  dependence  on 
the  parameters  1/a,  N(Ep),  and  lo,  while  Simanek  gives  no 
real  justification  for  the  values  he  uses.  Before  more 
definite  conclusions  can  be  reached  on  this  model,  accurate 
determinations  of  these  parameters  must  be  made. 
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6. 


CONCLUSION 


We  have  measured  the  extinction  by  small  particles 
o£  JVl  from  the  far  XR  to  the  far  UV.  This  was  done  because 
experimental  results  by  Granqvist  et.al.  (1976)  disagreed 
by  over  three  orders  of  magnitude  with  both  the  quantum 
mechanical  Gor'kov  Eliashberg  (1966)  theory  and  the 
classical  Drude  theory.  Our  measurements  served  to  confirm 
those  of  Granqvist  et.  al.  in  the  far  IR  and  extend  the 
data  to  the  UV.  Several  authors  had  proposed  theories  to 
account  for  Granqvist' s  observed  "anomalous"  absorption.  We 
have  extended  their  calculations  to  shorter  wavelengths  to 
compare  them  to  our  data.  Although  these  theories  show 
reasonable  agreement  in  the  far  IR,  they  differ  markedly 
from  our  experimental  data  at  shorter  wavelengths. 
Simanek's  (1980)  second  paper  is  a  possible  exception.  It 
shows  good  agreement  with  the  experimental  data  but  depends 
crucially  on  several  parameter.'*  which  are  not  accurately 
known. 

In  this  work  we  have  performed  a  shape  distribution 
calculation,  which  includes  the  effects  of  an  oxide  coating 
on  the  particles,  in  order  to  explain  the  experimental  daca. 
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This  calculation  is  based  on  the  classical  Drude  theory  and 
does  not  consider  quantum  size  effects.  The  shape 
distribution  in  our  model  is  a  uniform  continuous 

distribution  of  ellipsoids.  This  model  shows  excellent 
agreement  with  the  experimental  results  over  the  entire 
wavelength  region  measured.  In  order  to  determine  whether 
the  enhanced  absorption  is  due  primarily  to  shape  effects  or 
to  the  oxide  coating,  measurements  were  done  on  pure  aAl20j. 
In  the  far  IR  these  measurements  were  similar  to  those  of 
the  coated  Al  particles;  however,  for  wavelengths  shorter 
than  lOum,  the  absorption  by  the  pure  aAl20j  particles  is 
orders  of  magnitude  lower  than  for  coated  particles.  This 
would  imply  that  the  enhancement  is  due  to  shape  effects. 
On  the  other  hand,  it  is  possible  that  the  Mott  (1970)  a.c. 
conductivity  mechanism,  which  considers  the  effect  of 

defects  at  the  metal-oxide  interface,  enhances  the 
absorption  by  an  oxide  coating  as  compared  to  pure  aAl203 
particles.  This  is  the  basis  for  Simanek's  (1980)  second 
paper.  The  measurements  on  Al203  and  quartz  (Huffman  and 
Bohren  1980)  indicate  that  shape  effects  as  calculated  in 
the  CDE  model  are  sufficient  to  account  for  the  observed 
absorption.  Our  measurements  suggest  that  shape  effects  are 
also  the  major  contributing  factor  in  the  absorption  by 
oxide  coated  Al  particles,  but  further  measurements  should 
be  done  on  non-oxidizing  metals  to  ascertain  the  exact  role 
played  by  the  oxide  coating. 
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APPENDIX  A 


DETERMINATION  OF  DIELECTRIC  CONSTANTS 
FOR  POWDERED  MATERIALS 


In  Chap.  3  it  was  shown  that  the  shape  distribution 
model  (CDE)  agrees  well  with  the  experimental  results  for 
oA12Os,  thus  suggesting  the  validity  of  the  CDE  model. 
Assuming  this  model  works  equally  well  for  other  materials, 
we  have  developed  a  method  (developed  under  Naval  contract 
number  N00019-78-C-0479)  for  determining  IR  dielectric' 
constants  of  small  particle  systems  by  inverting  the  CDE 
calculation  with  the  aid  of  the  Lorentz  oscillator  theory. 
In  our  method  we  use  (29), 


a  "  Cahs/V  3  k  Im 


m 


along  with  the  expression  for  the  complex  dielectric 


constant  from  the  Lorentz  theory, 


_ 11  ,  .  2 

c 

Uj’-u’j-ll 


:a.  i) 


e^,  upij'  rj'  ai  1  wj  ace  parameters  in  this  method. 
N  is  the  number  of  oscillators  necessary  to  describe  the 
system?  usually  it  is  equal  to  the  number  of  peaks  in  the 
absorption  spectrum.  is  the  real  part  of  the  dielectric 
constant  at  infinite  energy,  and  normally  can  be  taken  to  be 
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the  square  of  the  handbook  value  for  the  index  of  refraction 
for  visible  light.  To  determine  the  dielectric  constants 
for  the  material,  the  fitting  parameters,  ok  ,  r^,  and  u>pij 
are  adjusted  via  an  interactive  computer  program  until  the 
calculated  absorption  ( acDE  )  is  in  good  agreement  with  the 
absorption  measured  experimentally  (otexp  ). 

This  technique  was  used  to  determine  the  dielectric 
constants  of  amorphous  AI2C)  (aAljOj).  For  amorphous 
materials  it  is  generally  not  obvious  what  value  of  N  to  use 
since  the  absorption  peaks  are  not  sharp  but  tend  to  be 
broad  and  unresolved.  In  this  case,  N  is  chosen  large 
enough  to  obtain  a  reasonable  fit.  For  aAl203  we  used  N»2. 
Fig.  A. 1  illustrates  the  fit  to  the  data  obtained  with  this 
method.  The  values  of  the  fitting  parameters  are  given 
under  the  figure.  The  dielectric  constants  are  found  using 
(A.l)  and  are  plotted  in  Fig.  A. 2. 

Our  method  of  determining  dielectric  constants  for 
small  particle  systems  is  very  similar  to  one  of  the 
standard  methods  for  the  determination  of  optical  constants 
from  reflectivity  measurements  on  bulk  solids.  For  this 
case,  Fresnel's  equation  is  used  rather  than  the  CDE 
equation  and  is  coupled  with  the  Lorentz  oscillator  theory. 
Although  the  dielectric  constants  obtained  by  our  method  may 
not  be  as  accurate  as  those  determined  from  reflectivity 
measurements  on  bulk  materials,  it  is  the  best  method 
currently  available  for  materials  which  are  only  available 

in  powder  form. 

Appendix  A 
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EXTINCTION  COEFFICIENT 


amorpn.  ALO 


FREQUENCY  (cm*1) 


Fig.  A-l  Extinction  Coefficient  of  aAl203 

Oscillator  fit  (open  circles)  compared 
to  measurements  (solid  line).  The 
oscillator  parameters  used  are; 

E°°*3’0' 

ui i *4 2 Ocm"* 1  ,  r, -150cm-1,  w*  -10  wf, 
to^-S^Ocm”* 1  ,  r  2*  160cm"  1  w*  *2.5  wj 


mmm 


FREQUENCY  (cm-') 


Dielectric  Constants  for  aAlj03 

Dielectric  constants  determined  using 
the  Lorentz  oscillator  theory  coupled 
with  the  CDE  theory 
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9 .  ABSTRACT 


Optical  properties  of  carbon  are  studied  in  bulk 
state  from  X  'v  0.05  to  100  um  for  graphite,  and  from 
X  'v-  0.05  to  1000  um  for  glassy  carbon;  in  small  particle 
state,  the  optical  studies  cover  the  spectral  range  going 
from  X  %  0 . 1  to  100  um  for  all  tho  materials. 

A  Krame rs-Kronig  analysis  of  near  normal  reflectance 
data  and/or  a  reflectance  data  fit  to  a  Dr ude-Lorentz  model 
gave  bulk  optical  constants.  Thesr  optical  constants  arc 
used  in  theoretical  calculations  of  extinction  and  the 
results  compared  with  experimental  results  obtained  from 
measurements  of  a  variety  of  carbon  particles.  It  is  in¬ 
ferred  that  the  high  experimentally  observed  extinction  is 
mainly  due  to  a  shape  effect. 


10.  INTRODUCTION 


Optical  properties  of  small  particles  have  become 
important  for  a  variety  of  reasons.  For  example,  atmospheric 
scientists  are  concerned  with  the  effects  of  small  particles 
on  weather.  It  is  known  that  the  role  of  atmospheric  aero¬ 
sols  in  the  earth's  heat  balance  depends  on  their  absorp¬ 
tivity  and  scattering  properties  (National  Academy  of 
Sciences,  1975).  In  order  to  predict  accurately  whether  the 
warming  tendency  due  to  absorption  and  the  cooling  tendency 
due  to  backscattering  in  the  earth's  atmosphere  will  dominate 
it  is  necessary  to  understand  the  optical  properties  of 
aerosol  particulates  (Twomey  and  Huffman,  1982)1. 

The  optical  properties  of  small  particles  are  impor¬ 
tant  in  astrophysics.  The  astrophysicist  needs  to  under¬ 
stand  the  extinction  by  small  particles  in  order  to  infer 
the  nature  of  the  solid  grains  responsible  for  the  observed 
interstellar  absorptions  (Huffman,  1977)  . 

Two  further  examples  demonstrate  the  great  diversity 
of  applications  of  small  particles.  For  military  applica¬ 
tions,  it  is  desirable  to  understand  the  optical  properties 
of  particles  wh*ch  are  to  be  used  to  generate  smoke  clouds 
for  obscuration  of  laser  weapons  and,  surveillance  devices; 
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for  certain  solar  energy  collector  types,  it  has  been  shown 
that  the  energy  collection  is  enhanced  by  doping  the  working 
fluid  with  very  small  absorbing  particles  (Hunt,  1980) . 

In  applications  of  the  later  kinds,  it  is  necessary 
to  maximize  the  absorption  by  small  particles.  Thus  carbon 
is  a  prime  candidate  material.  For  example,  Abdelrahman, 
Fumeaux,  and  Suter  (1979)  have  chosen  graphite  particles  as 
ideal  for  suspension  in  the  gas  of  a  solar  receiver.  Hunt 
(1978,  1979)  has  designed  a  collector-heat  exchanger  which 
will  utilize  carbon  particles  suspended  in  air. 

The  extinction  by  small  particles  can  be  calculated 
for  certain  shapes  if  the  optical  constants  as  a  function  of 
wavelength  are  known.  Although  the  literature  on  optical 
constants  of  carbon  is  abundant,  the  information  is  dis¬ 
parate  and  confined  to  limited  spectral  regions.  Workers 
rarely  agree  with  one  another  in  their  results.  (See  Twitty 
and  Weiman,  1971  for  a  review.) 

In  this  study,  two  carbonaceous  materials,  graphite 
and  glassy  carbon,  chosen  because  they  are  well  defined  (and 
thus  reproducible),  are  discussed.  Graphite  is  the  crystal 
form,  and  glassy  carbon  appears  to  be  the  best  defined 
carbon  variety  close  to  ideally  amorphous  carbon. 

Definitions  and  concepts  which  are  dealt  with  in 
this  work  are  covered  in  Chapter  11.  In  Chapter  12,  optical 
constants  of  graphite  and  glassy  carbon  are  surveyed: 
reflectance  data  of  these  materials  are  compiled  from  the 
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literature  and  augmented  by  our  own  measurements  and/or 
extrapolations.  The  data  are  analyzed  to  yield  optical 
constants  over  a  broad  spectral  range.  In  Chapter  13  ex 
tinction  by  fine  particles  and  smoke  of  carbon  is  experi 
mentally  measured,  and  results  compared  to  theoretical 
calculations  done  with  the  use  of  the  measured  optical 


constants . 


11.  REVIEW  OF  CONCEPTS  AND  DEFINITIONS 

< 

This  chapter  establishes  the  notation,  defines  the 
optical  constants  and  summarizes  the  theories  which  will  be 
used  in  this  work. 

After  defining  the  optical  constants  of  solids  and 
showing  their  relation  to  the  dielectric  constant,  the 
classical  dispersion  relations  and  the  Kramers-Kronig  rela¬ 
tions  are  reviewed.  The  extinction  by  small  particles  are 
introduced  through  the  Mie  theory  for  spheres  which  leads 
to  the  Rayleigh  approximations.  Rayleigh  theory  for  ellip¬ 
soids,  and  its  generalization  to  the  distributions  of 
ellipsoidal  shape,  are  then  discussed.  Finally,  a  sum  rule 
for  small  particle  extinction  is  introduced. 

11.1  Complex  Index  of  Refraction; 

Relation  to  the  Complex  Dielectric "Constant* 

Optical  properties  of  solids  deal  with  the  inter¬ 
action  of  electromagnetic  radiation  with  solid  natter.  This 
interaction  is  described  by  wavelength  dependent  complex 
index  of  refraction: 

m  ( >  )  =  n  *  i  k  ,  '  1  ' 


a  -let  ■«  i  led  treatment,  so*"*  any  E  f.  M  ‘ext,  Jack  sen 
:<  a  . 
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where  n  is  the  real  index  of  refraction  and  k  the  index  of 
absorption,  n  and  k  are  related  to  the  complex  dielectric 
function  e  through  the  relation: 


m2  ■  (n  +  ik) 2  ■  e 

From  e  “  Cj  +  ie2  it  follows  that: 

Ej  ■  n2  -  k1  (2) 

e2  ■  2nk  .  (3) 

The  (n,k)  or  ( e  x  /  c  2 )  pair  describes  completely  the  optical 
properties  of  a  material,  but  they  cannot  be  directly  mea¬ 
sured  by  experiment.  The  experimental  connection  is  made 
through  the  reflectance: 


R 


2 


where  r,  the  reflectivity  coefficient,  is  a  complex  function, 
defined  at  the  surface  of  a  solid,  as  the  ratio  of  the 
reflected  electric  field  Er  to  the  incident  electric  field 
.  At  normal  incidence: 


r 


(n-1)  +  ik 
(n+1)  +  i'K 


and 


R 


(n-1) 2  +  k2 

Tn+iyr;  "F" 


(4) 


(5) 
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11.2  Dispersion  Relations 


Classical  Dispersion  Relations* 

The  classical  theory  of  absorption  was  developed 
by  Lorentz  (for  insulators)  and  Drude  (for  free  electron 
metals).  The  combined  use  of  the  two  models  yields: 


C  M 


ID*  p  ID1, 

1 - y-lE, -  +  l  -j - £1-t - 

<DS  +  1YU!  .  j  Wj  “  ID*  —  lYjl" 


which  is  equivalent  to 


t  .  i  -pjH-"* 


^  “  UTtyT  +  j  (id  If i  +  Y  jiD^ 


+  l  ^£1 

A  t  m  * -7 

Ylfp . -v  +  I  "P3  w'j  ,  „ 

ID  (uTr+TJ!  j  (UjJ-ID3  )  !  +  Y  ?IDr 


«Y. 


(6) 


(7) 


where  ui,,  y j  and  wp^  are  respectively  the  resonance  frequency, 
the  damping  constant  and  the  plasma  frequency  of  the  jth 
oscillator.  They  determine  the  position,  the  width  and  the 
strength  of  the  jth  oscillator.  The  terms  out  of  the  summa¬ 
tion  signs  are  the  Drude  terms.  The  summation  is  over  the 
number  of  oscillators.  In  more  compact  form  (6)  and  (7) 
become : 


c 


i 


1  +  l 

j 


ID  2  ,  (  ID  *  -W  2  ) 

P  3  3 _ 

( ID  ?  “  U  2  )  2  +  YjWJ 


(8) 


*For  more  information,  see  for  example  Wooten  (1972,  Ch .  3). 


(9) 


y  “pfl  MYi  , 

4  fiJT-w*")  2  +  y  ^  w 2 


3 


The  Drude  terns  are  obtained  from  (8)  and  (9)  by  assuming 
that,  the  resonance  peak  of  the  first  oscillator  is  at  zero 


frequency  (w_,  ^  »  0)  . 


This  multiple  oscillator  model  is  used  in  the 
determination  of  optical  constants  of  solids.  From  a  reflec¬ 
tance  curve,  the  first  values  of  the  parameters  w.  ,  y^  and 


Up j  are  found  for  each  peak  by  locating  its  position,  and 


determining  its  width  and  its  strength.  They  are  used  in 
equations  (8)  and  (9)  to  calculate  e^tX)  and  e 2  ( X )  which, 
by  the  use  of  equations  (2)  and  (3)  yield  n(X)  and  k(X); 
the  n  and  k  obtained  are  used  in  (5)  to  yield  the  reflec¬ 
tance.  The  parameters  are  varied  until  the  calculated 
reflectance  fits  the  measured  reflectance  as  well  as  desired, 
giving  the  correspondincr  Cj(X)  and  e  2  ( X )  or  n  ( X )  and  k(X). 

This  technique  is  still  valid  even  if  the  solid  is 
anisotropic.  In  this  case,  the  sample  has  to  be  well 
oriented  and  polarized  light  must  be  used. 


The  Kramers-Kronig  Relations* 

The  reflectivity  coefficient  can  be  written  as 


r  (w) 


0  (  ui )  e 


16 


(10) 

(10) 


♦See  for  example  Wooten  (1972,  Ch .  6  and  Appendix  G) 
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.1. 1 


Equating  the  real  parts  and  the  imaginary  parts  of  (4)  and 
(6)  leads  to: 


n  *  1— *~-‘>2p-  co.  e  '  (U> 

,  2 p  s in  6  m 

k  *  r~p'*  •  2p  cos  5  ■  (12) 

The  knowledge  of  the  phase  0  will  allow  the  calculation  of 

t 

n  and  k.  The  determination  of  0  can  be  done  through  the  use 
of  the  Kramers-Kronig  relations,  which  connect  the  real  and 
imaginary  parts  of  the  response  function  of  a  linear  system, 
The  response  function  must  have  the  following  character¬ 
istics  : 

It  has  to  be  frequency  dependent. 

It  has  to  be  causal. 

It  has  to  be  analytic  in  the  upper  half  of  the 
complex  u>  plane. 

It  can  be  shown  that  the  reflectivity  coefficient  is  a 
response  function  between  the  incident  and  reflected  wave 
at  the  surface  of  a  solid.  It  can  also  be  shown  that 

In  r  (id)  *  In  p  (id)  +  16  (id) 

satisfies  the  requirements  for  the  application  of  the 
Kramers-Kronig  relations. 

The  relation  between  0  (id)  and  c(u)  which  is  cf 
interest  in  this  treatment  is: 
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s  <w) 


'o 


In  o  (ui 1 ) 

Ul  !  -_0)2 


dw 1 


where  P  is  a  principal  part.  The  removal  of  the  principal 
part  and  the  replacement  of  p(ui’)  by  |R(w')  | ^  gives: 


e  (u) 


W  r  m  [R(g>M/R(a»)3  duJ. 

TT  0  U2  ~  Ui '  2 


In  terms  of  wavelength,  the  substitution  of  us  »  2irc/X  can  be 
made  to  yield: 


9  (X) 


2X 


m  [»»■)/»[»]  ^ 


i  _ 


(13) 


Relation  (13)  suggests  that  the  reflectance  has  to 
be  measured  from  2ero  to  infinite  wavelength.  This  is 
impossible  in  practice;  in  general,  the  reflectance  is 
measured  for  X  such  that: 


X  .  *  X  £  X 

mm  max 


Thus,  (13)  can  be  rewritten: 


e  (X) 


2X 

TT 


(min  In  [R ( X • ) /R ( X ) 1  , 

— -  d* 


+  2  X  fAmax  In  [  R  (  X  1  )/R(X)  ] 

r  ;>  x '  ~  x2  ' 

min 


dX'  +  ^ 


\ 


In  [r ( >  1  ) /R  ( X ) J  . 

’  i :  t  a 


T' 


max 


( x )  +  e .  ( ■ )  +  e ,  ( x ) 


e2  is  calculated  with  the  actual  experimental  data.  For 
OjOO  and  9, (A),  some  sort  of  extrapolation  has  to  be  done. 
For  metal  and  metal-like  materials,  the  Drude  model  can  be 
used  to  £ind  R(A')  and  compute  9,  (A).  For  other  kinds  of 
materials,  Steyer  (1974)  derived  a  relation  based  on  the 
classical  dispersion  theory. 

The  short  wavelength  spectral  range  is  more  diffi¬ 
cult  to  treat.  To  include  all  the  interband  transitions, 
the  reflectance  is  approximated  by; 

R  ( A )  -  R.  x  ( A/A  )S  ,  (14) 

le  le 

where  A  <  A,  and  s  >  0.  R,  and  A.  are  the  reflectance 
le  le  le 

and  the  wavelength  at  the  last  experimental  point.  At  lower 
wavelengths,  the  reflectance  is  approximated  by: 

R(A)  -  Rf  x  (A/>.£ )  "  .  (15) 

Rj  is  calculated  by  the  use  of  (14)  at  the  point  where  the 
second  extrapolation  starts. 

With  this  extrapolation,  it  is  assumed  that  the 
wavelength  is  low  enough  to  allow  a  f ree-electron  descrip¬ 
tion  of  the  electrons. 

These  extrapolations  are  used  in  the  integration  of 
9 . ( \ )  to  yield  (Wooten,  1971,  D.  249): 


e,  (X) 


l 

I? 


In 


,  f 

imr. 


■]  -  M 


Jr 


r  (Xle/X)2n+1 
“  T5n+TP 


le* 


+  (4-s) 


(X£/X) 


2n+l 


2n+l) 2 


]• 


(16) 


The  parameter  s  and  the  wavelength  X£  are  chosen  such  that 
the  calculated  optical  constants  are  in  agreement  with  those 
obtained  from  an  independent  measurement  in  a  limited  spec¬ 
tral  range. 

The  actual  reflectance  .data  and  the  extrapolations 
are. used  in  the  numerical  integration  which  is  done  by  the 
use  of  Simpson's  rule  of  integration. 


11.3  Extinction  by  Small  Particles . 

Introduction  to  Mie  Theory 

A  beam  of  light  which  passes  through  a  collection  of 
particles  is  attenuated.  This  attenuation,  due  to  scatter¬ 
ing  and  absorption  by  the  particles,  is  called  extinction 
which  is  defined  to  be: 


Extinction  *  Scattering  +  Absorption  , 

A  way  for  inferring  the  extinction  is  to  calculate 
the  transmission  T  of  the  beam  through  a  path  length  £  of 
the  distribution  of  particles.  For  N  spherical  particles 
per  unit  volume,  each  of  radius  a: 
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T 


exp  (-NQextTra2 )  l 


(17) 


I/I 


where : 


®ext  "  ^sca  +  ®abs 


The  Q's  are  called  efficiency  factors  for  spheres,  and  are 
related  to  the  following  quantities: 

Csca  "  ^sca  *  a*  is  the  cr°8S  section  for  removal 

of  energy  from  the  incident  beam  by  scattering, 
cabs  "  ^abs  *  ft2  tlle  cr0SB  section  for  removal 

of  energy  from  the  incident  beam  by  absorption, 


C 


ext 


c 

sea 


+  c 


abs 


(18) 


Mie  (1908)  and  Debye  (1909)  developed  a  theory  for  calcu¬ 
lating  the  efficiency  factors  for  spheres  of  arbitrary 
sizes.  This  theory,  usually  referred  to  as  the  Mie  theory 
is  presented  in  books  by  van  de  Hulst  (1957),.  Kerker  (1969) 
and  Bohren  and  Huffman  (1983). 

From  the  theory  s 

Q  -  2  l  (2n+l)  [  |a  |  2  +  |b  ! 2]  ,  (19) 

sea  x*  n«i  n  n 

Qabs  *  7*  I  (2n+1)  Re(an  +  bn-  '  (20) 

nrel 

where  x  «•’  2ua/X  is  the  size  parameter.  The  coefficients 
an  and  bn  are  complex  expressions  involving  spherical  Bessel 
and  Hankel  functions  and  their  derivatives.  They  can  be 
computed  for  given  values  of  the  complex  index  of  refraction 
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n  and  the  size  parameter  x.  The  computation  procedure  is 
given  by  Wichramasinghe  (1973)  and  related  computer  programs 
are  in  the  book  by  Bohren  and  Huffman  (1983) . 

Theory  and  experiment  are  connected  as  follows. 

If  we  define: 


N  «  n/V 


where  n  is  the  total  number  of  particles  in  the  volume  v, 

C  ■  m/V  *  ^  ^  tt  a 3  p  , 


where  m  is  the  total  mass  of  particles  and  d  the  mass  den¬ 
sity  of  the  particle  material,  (17)  becomes: 

I/Io  -  exp  (-  a 

«=  exp  (-  a  |)  ,  (21a) 

a  -  S.  In  (I0/I)  .  (21b) 

If  the  optical  density  (O.D.)  defined  as:  O.D.  *  277  ln  ~y 
is  used,  then  a  ■  2 . 3  ~  (O.D.),  where  c  is  the  mass/area  of 
particles.  Also 


0  .  ti  a 

ext 

4/3  if  a1" 


3  ^ext 


a  is  the  volume  normalized  extinction  coefficient  and  is 
related  to  the  mass  extinction  coefficient  a  by: 


In  the  case  of  a  size  distribution  of  spheres: 


( 


N  (a)  Q(a)  rra3da 

a  =  1 - -  -  - ...  t 

N  (a)  ^  Tra3  da 

where  N(a)  is  the  number  of  particles  per  unit  volume  with 
radii  varying  between  a  and  a  +  da. 

From  (21) ,  a  can  be  determined  experimentally  by 
measuring  the  transmission  and  c.  The  use  of  (22)  leads  to 
an  experimental  value  of  Qext- 

If  the  particles  are  very  small  compared  to  the 
wavelength  (a  <<  A )  and  when  the  phase  shift  of  light  in  the 
particle  is  negligible  (|m|x  <<  1)  (Rayleigh  approximation), 
then  the  Mie  efficiency  factors  can  be  expanded  in  power 
series  as  functions  of  the  si2e  parameter  x;  the  lowest 
power  terms  of  importance  in  a  first  approximation  are: 


0  ® 

8 

3  * 

h2"1! 

2 

(24) 

V 

sea 

vm2+?, 

/ 

^abs 

4  x  Im 

\m2  +  ^ 

} 

• 

(25) 

Relations  (24)  and  (25) 

apply  when 

the 

particles 

are  in  a 

medium  with  complex  index  of 

refraction 

m' 2  =  1. 

They  can 

be  generalized  to  give: 

^sca 

8  x1*  1 

m 

2 

(26) 

3  X  1 

g  +  2c 

nr 

t 

Qabs 

r€-Z 

a  v  Tm  J  m 

A 

(27) 
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where  c  is  the  dielectric  constant  of  the  particles  and  e 

TTl 

is  the  dielectric  constant  of  the  medium  in  which  they  are 
embedded. 

The  Rayleigh  approximations  have  introduced  two  main 
simplifications  in  the  treatment  of  extinction  by  small 
particles.  First  of  all,  there  is  no  shape  effect  in  volume 
normalized  extinction,  as  illustrated  in  Fig.  17  •  At 
small  radii,  the  extinction  is  a  constant.  However,  for 
radii  above  0.012  ym,  the  extinction  is  extremely  radius 
dependent.  One  notices  a  decrease  in  extinction  for  radius 
larger  than  0.1  ym.  This  decrease  is  due  to  saturation 
effect?  when  absorption  is  very  effective,  the  inner  parts 
of  the  sphere  do  not  participate  in  the  absorption  process. 
With  larger  spheres,  more  of  the  inner  volume  is  ineffective, 
which  results  in  a  decrease  of  the  extinction  per  unit 
volume. 

From  an  experimental  point  of  view,  one  does  not 
have  to  worry  about  the  size  distribution  of  particles  if 
they  are  small  enough  compared  to  the  wavelength.  This 
occurs  easily  in  the  IR  spectral  range. 

The  second  simplification  introduced  by  tne  Rayleigh 
approximations  is  that  the  shape  effects  can  be  more  easily 
treated  theoretically,  as  will  be  shown  in  the  next  section. 


Radius  (^m) 


Absorption  and  scattering  cross  sections  normal 
ized  per  unit  volume  of  solid  for  carbon  in  the 
visible,  calculated  using  Mie  theory 


11  * 4  Rayleigh  Theory  for  Ellipsoids;  The  Shape  Distribution. 


Small  particles  have  been  discussed  with  the  assump¬ 
tion  that  they  are  spherical.  Although  this  is  not  very 
often  true,  the  spherical  hypothesis  allowed  the  use  of  the 
Mie  theory  which  can  be  generalized  to  infinite  cylinders 
only.  A  more  realistic  shape  for  particles  is  the  ellipsoid, 
for  which  polarizability  in  an  electric  field  can  be  deter¬ 
mined  and  used  in  the  derivation  of  an  approximate  expres¬ 
sion  for  the  absorption  and  scattering  cross-sections,  if 
the  Rayleigh  conditions  are  valid.  The  derivation  of  the 
polarizability  can  be  found  in  the  bibliography  given  for 
the  Mie  theory. 

For  an  ellipsoidal  particle  with  semi-axes  a,  b  and 
c  in  an  electric  field  parallel  to  the  axis  (j  *  1,2,3)  , 
polarizability  is  given  by: 

aj  "  TF  e'-'+L-TF-ry  '  (28) 

J  m  3  m 

where  V  =  (4/3  )ir  abc  is  the  volume  of  the  ellipsoid.  The 
depolarization  factors  Lj  are  given  by: 

L  -  S§£  f - ^ c  ,  (29) 

2  J  (s+a2)  { (a2+s) (b2+s) (c2+s) P 

where : 


s  =  a  +  b  +  c 
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L2  and  I,3  are  computed  the  same  way,  with  the  right  substi¬ 
tution  in  (29).  Using  the  integrals  by  which  the.  L^'s  are 


calculated,  it  can  be  shown  that: 


3 

I  L 
j-1 


'J 


-  1 


(30) 


In  the  Rayleigh  approximations,  Csca  and  C^g  for  a  single 


ellipsoid  of  abitrary  Lj  ares 


k- 


e-e 


sea 
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'abs 


kV  Im 


(31) 


(32) 


The  superscript  j  shows  that  E  field  is  parallel  to  the  jth 
axis;  k  «  2-tt/ X .  For  a  collection  of  ellipsoids  randomly 
oriented,  an  average  value  is  taken;  the  absorption  case 
is  given: 


<C  K  > 

abs 


l  k 


e  t 
m  2 


3  j.!  Lj  1  le  j  +  c^d/Lj-l)  J  2  +  z\) 


33) 


From  (33)  it  can  be  seen  that  resonance  occurs  whens 


Li£i  * 


or 


- 


:  34 ) 
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Since  0  <  L  <  1,  the  resonance  occurs  only  for  negative 

j 

values  of  e  t .  For  spheres,  Lj  »  1/3  for  all  directions,  and 
the  resonance  is  obtained  for: 


e,  *  -2e 

1  in 

A  distribution  of  identical  ellipsoids  with  arbitrary  Lj ’ s 
will  give  three  absorption  peaks  because  of  the  three  values 


of  L j . 


A  more  general  treatment  of  the  preceding 


case  is  done  by  considering  a  variety  of  shapes  of 
ellipsoids . 

Given  a  shape  distribution  expression  P(Lj,Lz),  it 
is  possible  to  compute: 


<<C  .  )> 

abs 


?U,,LZ)  <C,b  >  dL,  dL; 


The  integration  is  done  over  Lj  and  L2  only  since  from  (30) 
it  is  seen  that  only  two  L ^ 1 s  are  independent. 

Assuming  that  all  shapes  of  ellipsoids  are  equally 
probable,  which  mathematically  means  that  P(Li,Lj)  is  a 
constant,  Huffman  and  Bohren  (1980)  found  (35)  to  be: 


<  <  C  L  >  >  -  k<  V)  Im 


log  (c/e, 


.)!  ■ 


(For  the  log  calculation,  use  lej  and  — ir  <  0  <  n . )  For  this 
treatment,  a  system  of  particles  with  all  shapes  of  ellip¬ 
soids  should  produce  absorption  over  the  who '  .  spectral 
range  where  e1  is  negative.  Rathmann (1981 ,  p.  18)  calcu¬ 
lated  aluminum  absorption  coefficients  for  spheres,  single 
ellipsoids  and  a  continuous  distribution  of  ellipsoids. 

Relation  (36)  is  used  in  our  theoretical  calcula¬ 
tions,  and  the  corresponding  curves  are  labeled  CDE 
(continuous  distribution  of  ellipsoids) . 

Sum  Rule  for  Extinction 

Bohren  and  Huffman  (1983)  have  derived  a  sum  rule 
for  extinction  by  spherical  particles: 

WJ>  -  4»Js'  ejjSC-iVS  •  w 

For  conducting  materials,  e(X"«)  is  large;  thus  (37)  has  an 
upper  limit  such  that: 

,°° 

CextU)  dX  <  4tt  5a5  .  (38) 

o 

Using  (22)  in  conjunction  with  Cext  *  gives: 

Cext  (X)  ■  ^  ira3  a(X)  .  (39) 
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A  substitution  of  (39)  in  (38)  leads  to: 


/<*> 

o(X)  dX  1  3tt2 
o 

This  relation  will  be  used  to  check  the  sphericity  of 
experimentally  produced  particles. 

The  sum  rule  (37)  was  first  derived  by  Purcell 
(1969)  for  spheroidal  particles. 


12.  BULK  OPTICAL  CONSTANTS 


Carbon  is  the  generic  designation  of  a  large  variety 
of  carbonaceous  materials  including  graphite  and  amorphous 
carbon.  Graphite  and  amorphous  carbon  are  well  defined,  but 
many  types  of  coal  and  soot  differ  in  composition  and  struc¬ 
ture.  According  to  Dalzell  and  Sarofin  (1969) ,  soots  are 
randomly  mixed  graphitic  particles  in  a  matrix  of  amorphous- 
carbon;  they  differ  from  one  another  by  their  hydrogen/ 
carbon  ratio.  The  actual  difference  in  the  composition  of 
carbonaceous  materials,  and  apparent,  differences  resulting 
from  the  use  of  inappropriate  experimental  technique  in  the 
determination  of  the  optical  constants,  have  led  to  results, 
which  in  some  cases,  are  hardly  comparable.  Disagreements 
in  results  are  shown  through  the  review  of  some  previous 
work  presented  below.  The  large  variation  in  published 
optical  constants  can  be  seen  in  Pig.  18. 

Stull  and  Plass  (1960)  derived  dispersion  relations 
for  carbon,  which  they  used  to  fit  some  other  workers' 
reflectance  data  between  \  ^  0.436  urn  and  13  urn.  The  same 
method  has  been  used  by  Dalzell  and  Sarofin  (1969)  to  fit 
soot  pellet  reflectance  data  in  about  the  same  spectral 

104 


-•‘—Arc Vaporized  Carbon(2) 

- Arc  Vaporized  Carbon  (1 ) 

- Glassy  Carbon 


9  E 


^  C  vo 
cn  m  o\ 
r '  H  H 
<j\  0 
H  U  - 
£ 

-•O  4J 

ns  c  u 
j  ns  b 

io  t 
X  0  0 

to  tj  as 

M  M 
<  10  tJ 
2  C 

<o  **h  n 
c  Q 

IQ  u 
T3  0) 
W  C  4-1 

E  ro  in 
io  o 

•iH  In 


c 

H  H  w 
iH  ' — 

in 

*H  iH 

4J 

i  ^  «i 

c 

—  0  0 

IQ 

Hh  U 

4J 

e 

w 

0t',O 

c 

£  r»  C 

0 

V4  OV  Ifl 

0 

IQ  r-4 

0  0) 

H 

»  4J 

(0 

>i  1  'H 

0 

in  H  £ 

•H 

in  (Q  tt. 

4J 

|Q  (0 

Qj 

pH  4J  l-l 

O 

O  0>  O* 

'ks 

is 

i>2 

[»%> 

N-«l 

t2& 

ii 

LiJ 

range.  The  results  of  the  two  groups  do  not  agree  with  one 
another.  More  recently,  Lee  and  Tien  (1981)  derived  a  dis¬ 
persion  model  based  on  electronic  band  structure  considera¬ 
tions.  The  dispersion  constants  were  determined  from  the 
transmission  data  of  soot  flame.  They  questioned  their  pre¬ 
decessors'  methods  and  results.  Tomaselli  et  al.  (1981) 
measured  reflectance  of  pressed  pellets  of  various  carbona¬ 
ceous  materials.  They  graphically  analyzed  data  by  plotting 
iso-reflectance  curves  and  reading  n  and  k  from  published 
tables.  None  of  their  results  came  close  to  published 
results  for  similar  materials.  DiNardo  and  Goland  (1971) 
derived  optical  constants  of  a,rc  evaporated  carbon  film  by 
matching  the  transition  radiation  theory  to  experiments  from 
X  ^  0.23  urn  to  0.56  pm.  Arakawa ,  Williams,  and  Inagaki 
(1977)  studied  a  similar  material  from  X  t  0.33  um  to  2.1 
um..  A  comparison  of  the  two  results  shows  a  clear  disagree¬ 
ment  (see  Pig. 18.  ).  According  to  Arakawa  et  al.  (1977), 
the  optical  properties  of  arc  evaporated  carbon  vary  with 
the  conditions  of  preparation  of  the  film.  Pluchino  et  al. 
(1980)  isolated  a  single  micro-sized  particle  of  carbon; 
they  electrostatically  suspended  and  irradiated  it  with  a 
laser  beam.  The  scattered  radiation  intensity  was  measured 
as  a  function  of  angle.  Electron-microscopy  showed  that  the 
particle  was  spherical.  Data  were  then  analyzed  with  the 
Mie  theory.  At  X  'v  0.48  um  and  for  a  particle  radius  a  = 
3.06  um,  n  and  k  are  found  to  be  1.7  and  0.8  respectively. 
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This  review  along  with  Fig.  18  has  shown  that 
there  is  no  general  consensus  on  optical  constants  of 


carbon;  it  has  shown  also  that  available  optical  constants 
are  confined  to  a  limited  spectral  range. 

In  this  section,  we  survey  the  optical  constants  of 
graphite  and  glassy  carbon,  and  gather  existing  reflectance 
data  on  these  materials.  In  order  to  calculate  the  optical 
constants  over  a 'wide  spectral  range,  the  reflectance 
information  gathered  is  completed  by  our  measurements  and/or 
extrapolations  wherever  it  is  necessary,  to  yield  a  reflec¬ 
tance  curve  over  a  wide  spectral  range ,  for  use  in  the 
Kramers-Kronig  analysis.  A  fit  of  the  reflectance  data, 
based  on  the  combined  use  of  Drude  and  Lorentz  dispersion 
relations  has  also  been  done. 

Techniques  for  the  determination  of  optical  con¬ 
stants  are  given;  then  graphite  and  glassy  carbon  are 
studied  in  turn. 

12.1  Techniques  for  Measuring  Optical  Constants* 

Optical  properties  of  solids  are  due  to  the  wave¬ 
length  dependence  of  n  and  k  (£j  and  e2)  which  cannot  be 
measured  directly  experimentally.  The  techniques  used  to 
collect  information  leading  to  optical  constants  are  sum¬ 
marized  in  Table  3. 


Table  3.  Techniques  for  determining  optical  constants  of 
solids 

Prom  Huffman  (1977) 


Technique 

Cc  ents 

Transmission  methods 

(1) 

n  from  minimum  deviation? 
k  from  transmission  and  n 

High  accuracy  for  n 

(2) 

n  and  k  from  transmission 
and  normal  incidence  reflec¬ 
tance  measurements 

Easy  but  useful  only  in 
relatively  transparent 
regions 

Reflection  methods 

(1) 

Two  polarized  reflectance 
measurements  at  one  oblique 
incidence  angle 

Larger  sample  sizes 
required 

(2) 

Two  reflectance  measurements 
at  different  angles 

(3) 

Determination  of  special 
angles  (i.e.,  Brewster's  angle! 
and  reflectance  there 

) 

(4) 

Determination  of  ellipso- 
metric  parameters 

Sensitive  to  surface 
films 

(5) 

Reflectance  in  a  broad 
wavelength  range  and 
Kramers-Kronig  analysis 

Extrapolation  to  0  and 
to  00  necessary 

(6) 

Reflectance  vs.  wavelength 
and  oscillator  model  fit 
to  the  data 

Compact  presentation  of 
results  as  oscillator 
parameters 

Other  methods 


(1)  Electron  energy  loss 
measurements 


Does  not  require 
polished  surfaces 


12.2  Graphite . 


Perfect  graphite  crystals  are  rare.  Figure  19. 
represents  the  ideal  natural  graphitic  structure.  The  basal 
planes  are  such  that  every  other  plane  fits  exactly  over  the 
first;  this  is  the  ABAB  arrangement.  Another  one  is  such 
that  the  basal  planes  are  in  the  ABC ABC  sequence. 

The  atoms  of  carbon  are  associated  in  hexagonal 
rings.  Each  atom  has  four  electrons  on  its  outer  shell; 
three  of  them  (o-elect.rons)  participate  in  holding  together 
the  hexagonal  ring  by  covalent  bonds;  the  fourth  one  (tt- 
electron)  moves  in  the  ring  and  contributes  to  the  bonding 
between  planes.  Many  properties  of  graphite,  including  its 
optical  properties,  are  highly  anisotropic.  The  anisotropy 
of  the  optical  properties  is  due  to  the  ease  with  which  the 
ir-electrons  move  in  the  layer  planes;  higher  energy  is 
needed  to  jump  from  plane  to  plane. 

Forms  most  closely  approaching  the  ideal  structure 
of  the  material  graphite  are  pyrolytic  graphite  and  annealed 
pyrolytic  graphite.  Pyrolytic  graphite  is  produced  by  a 
decomposition  of  a  hydrocarbon  on  a  hot  surface.  Although 
it  is  not  a  true  graphite,  its  properties  are  extremely 
anisotropic.  It  can  be  converted  into  almost  perfect  single 
crystals  by  heating  at  temperatures  above  2900 °C  for  several 
hours.  It  is  then  called  (stress)  annealed  pyrolytic 
graphite  (Shobert,  19  64)  . 


The  optical  properties  of  graphite  depend  on  the 
direction  of  the  E-field  with  respect  to  its  optical  c-axis, 
E  JL  C  and  E  U  C  leading  to  and  e/;  respectively.  Because 
samples  can  easily  be  prepared  by  cleavage  of  the  crystal 
in  the  basal  plane ,  many  optical  constants  for  the  E  i  C 
direction  are  available  in  the  literature.  Fewer  measure¬ 
ments  have  been  done  for  E  //  C  since  the  sample  preparation 
is  considerably  more  difficult  in  this  case.  El  C  and 
E  ft  C  optical  constants  of  graphite  will  be  reviewed  in 
turn. 

Survey  of  Optical  Constants  of  Graphite 

Many  workers  measured  reflectance  on  natural  graph¬ 
ite,  pyrolytic  graphite  or  annealed  pyrolytic  graphite. 

They  analyzed  their  data  by  the  Kramers-Kronig  method  (Taft 
and  Philipp,  1965;  Philipp,  1977),  by  fitting  (Greenaway  et 
al.,  1969;  Klucker,  Skibowski,  and  Steinman,  1974)  or  by  the 
use  of  Fresnel's  equations  (Carter  et  al. ,  1965). 

Zeppenfeld  (1967)  was  the  first  to  use  the  electron 
energy  loss  method  for  deriving  optical  constants  of  graph¬ 
ite.  This  method  consists  of  determining  the  properties  of 
bulk  materials  by  measuring  the  transmission  of  electron 
beams  through  thin  films,  t x  and  e2  are  obtained  by  the 
Kramers-Kronig  analysis  of  the  energy  loss  function  f(u>) 
which  is  related  to  the  dielectric  function  by  (Daniels  et 
al. ,  1970) : 
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f  (w)  *  -  Im  (1/e) 

The  survey  of  optical  constants  is  summarized  in  Table  4 . 

Comparison  of  results  for  E  J.  C  is  seen  in  Fig.  20. a 
which  shows  e2  in  a  spectral  range  which  contains  a  strong 
peak.  Good  agreement  exists  in  the  shape  of  the  curve,  but 
the  discrepancy  in  the  magnitude  of  the  peak  is  obvious. 

For  E  //  C,  Fig.  20. b  shows  complete  disagreement 
between  the  results  of  the  reflectance  techniques  and  those 
of  the  electron  energy-loss  method.  Because  of  the  severity 
of  the  discrepancy,  Venghaus  (1975)  repeated  the  electron 
energy-loss  measurements  on  graphite.  His  results  were  in 

agreement  with  the  measurements  ..done  with  the  same  technique 

(2  ,  T  1 

by  other  workers^  but  still  in  disagreement  with  reflectance 
method  results.  Wessjohan  (Venghaus,  1975)  computed  reflec¬ 
tance  using  Venghaus  (1975)  and  Klucker  et  al.  (1974) 
e  data;  Venghaus 's  results  compared  favorably  with  experi¬ 
ments,  but  those  from  Klucker  et  al.  failed  to  do  so.  Al¬ 
though  Klucker  (1971),  who  noticed  the  discrepancy  between 
his  experimental  results  and  calculations,  assigned  this 
failure  to  scattered  light,  Venghaus  concluded  that  the 
electron  energy-loss  method  gives  better  results. 

In  IR  (X  %  2.5  um  to  100  urn),  Venghaus  measured 
the  reflectance  on  graphite  in  the  E//C  direction  and  ana¬ 
lyzed  the  data  by  the  Kramers-Kronig  method. 
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Survey  of  optical  properties  of  graphite 


-Continued 


As  computed  by  Klucker  et  al.  (1974) ,  Taft  and 
Philipps  (1965) ,  and  Zeppenfeld  (1967) 


This  Study's  Investigation  on  Graphite 

In  order  to  have  a  coherent  set  of  optical  constants 
in  a  broad  spectral  range  for  the  E  i  C  direction  of  graph¬ 
ite,  and  also  to  investigate  the  discrepancy  of  the  c2  peak 
at  A  0.088  pm,  Philipp's  (1977)  extended  reflectance  data 
(see  Fig.  21.)  have  been  analyzed  by  the  Kramers-Kronig 
method.  Philipp  obtained  his  extended  reflectance  data  by 
combining  previously  published  results  (Taft  and  Philipp, 
1965),  measurements  between  A  ^  12,4  urn  and  41  pm,  and 
Sato's  (1968)  data. 

It  is  known  that  the  Kramers-Kronig  analysis  results 
depend  strongly  on  the  extrapolation  at  low  wavelengths 
(Wooten,  1972).  To  satisfy  the  requirements  of  the  extrapo¬ 
lation,  which  is  explained  in  Chapter  11,  some  guide  values 
of  n  and  k  found  from  an  independent  experiment  are  needed. 
Huffman  (1979) ,  in  work  done  at  the  University  of  Wisconsin 
monochromatic  radiation  facility,  measured  the  transmission 
of  a  0.27  pm  thick  cleaved  graphite  sample  from  A  1.25  pm 
to  0.15  pm.  The  values  of  n  and  k  found  from  this  experi¬ 
ment  have  been  used  as  tests  for  the  choice  of  the  parameter 
s  and  the  wavelength  A.  These  values  are  s  ■  2,  A ^  %  0.0258 
pm.  Taft  and  Philipp  used  the  results  of  the  transmission 
of  a  400A  thick  sample  from  A  v  0.113  pm  to  0,177  pm  as 
guiding  values  in  their  calculations.  Results  of  our  calcu¬ 
lations  are  shown  .in  Figs.  22.,  23  .,  and  24.  and  in 
Appendix  B. 
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Fig.  22.  Real  part  of  the  El  C  polarization  of  graphite 
as  computed  by  the  Kramers-Kronig  method 
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Fig.  24 .  Optical  constants  n  and  k  for  the  El  C  polarization  of  graphite 
computed  by  the  Kramers-Kronig  method 


The  reflectance  data  fit  result  is  shown  in  Fig.  25. 

A  discussion  of  the  discrepancy  observed  between  the  data 
and  the  fit  from  X  =  0.3  to  15  pm  is  given  in  the  "Discus¬ 
sion  of  Results"  section  of  this  chapter.  The  oscillator 
parameters  used  are  presented  in  Table  5. 

Representative  values  of  e  t  and  e2  for  the  E //  C 
direction  of  graphite  derived  from  the  reflectance  and  the 
electron  energy- loss  measurements  have  been  converted  to  n 
and  k  and  plotted  (Fig.  26.  )  .  Venghaus's  results  below 
X  v  2,5  m  were  obtained  privately.  Venghaus  (1977)  argued 
that  reflectance  measurements  are  not  seriously  affected  by 
the  roughness  of  the  reflectance  surface  in  the  IR  wave¬ 
length  range;  n  and  k  obtained  by  optical  methods  are  thus 
reliable  from  2,5  pm  to  100  pm.  In  the  short  wavelength 
range  (X  v  0.03  to  about  0.22  pm),  the  increasing  dependence 
on  the  reflectance  on  the  state  of  the  surface  as  the  wave¬ 
length  decreases  makes  the  results  of  the  optical  method 
less  reliable  than  those  obtained  with  the  electron  energy- 
loss  technique.  Since  above  X  v  0.25  pm,  the  Cerenkov 
radiation  gives  the  dominant  contribution  to  the  electron 
energy-loss  for  E  //  C  (Tossati  and  Bassani,  1970),  the  elec¬ 
tron  energy-loss  results  have  been  disregarded  from  X  v 
0,25  pm  to  2.5  pm.  From  0.15  pm  to  0.25  pm,  geometrical 
considerations  have  been  used.  The  slope  of  the  electron 
energy-loss  n  curve  at  short  X  ^  0.15  um  is  the  same  as  the 
slope  of  the  reflectance  n  curve  from  X  'u  0.21  pm  to  0.2  5  um. 
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Table  5,  Oscillator  parameters  for  the  fit  of  the 
reflectance  curve  for  the  £j_  C  direction 
of  graphite 


Parameter 

Value  (cm-1) 

Oscillator  1 

WP1 

1200 

W1 

0.00 

Y1 

15.00 

Oscillator  2 

WP2 

1300 

W2 

120 

Y2 

120 

Oscillator  3 

WP  3 

6  x  103 

“3 

4  x  103 

Y3 

7.5  x  103 

Oscillator  4 

WP4 

110  x  10 3 

W4 

98  x  103 

Y4 

25  x  103 

123 


ErQun  Graphite.  E^C 
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Optical  constants  n  and  k  of  the  E  //  C  polarization  of  graphite 


Based  on  that  consideration,  a  connecting  curve  has  been 
drawn  between  the  two  n  curves  from  X  ^  0.15  pm  to  0.21  pm. 
Prom  this  combination  of  optical  method  results  and  elec¬ 
tron  energy-loss  measurements,  a  set  of  n  and  k  data  for 
E //  C  has  been  put  together  from  X  ^  0.025  pm  to  100  pm. 

(See  Figs.  27.  and  28..  and  Appendix  B.  ) 

For  the  E //  C  direction  of  graphite,  the  reflectance 
computed  from  n  and  k  extracted  from  the  literature  has  been 
fitted  using  a  Dr ude-Lo rents  oscillator  model;  the  result 
is  shown  in  Fig.  29.'  The  oscillator  parameters  used  are 
presented  in  Table  6. 

12.3  Glassy  Carbon . 

Glassy  carbon,  also  known  as  vitreous  carbon,  is 
produced  by  the  slow  pyrolysis  of  polymers  such  as  cellu3en 
and  aromatic  resins.  It  looks  like  black  glass,  has  a  low 
density  (.1.4 -1.5  g/cm1)  and  also  a  low  porosity  and  perme¬ 
ability  to  water.  It  has  about  the  same  hardness  as  glass; 
it  withstands  high  temperature  in  absence  of  oxygen.  X-ray 
studies  have  shown  that  glassy  carbon  has  randomly  oriented 
crystals  which  are  not  very  much  altered  by  heat  treatment. 
At  3273°K,  glassy  carbon  crystallite  size  is  6  mm?  they  are 
smaller  at  lower  temperatures  (Halpin  and  Jenkins,  1969). 

Taft  and  Philipp  (1965)  published  near  normal  inci¬ 
dence  reflectance  data  of  glassy  carbon  up  to  26  eV. 

Williams  and  Arakawa  (1972)  published  plots  of  n  and  k 


Adjusted  optical  constants  n  and  k  of  the  E  //  C  polarization  of  graphite 


Fig  -  29.  Reflectance  of  the  E  if  Q.  polarization  of  graphite 


Table  6. 


Oscillator  parameters  for  the  fit  of  the  reflec¬ 
tance  curve  for  the  E  //  C  direction  of  graphite 


Oscillator  1 


Parameter 


Value  (cm”1) 

8500 

0.00 

900 


Oscillator  2 

WP2 

2950 

U2 

1500 

Y2 

1500 

Oscillator  3 

WP3 

5,000 

W3 

25,000 

Y3 

3,000 

Oscillator  4 


w 


P4 


95,000 


i 


of  glassy  carbon  from  0  to  82  eV,  but  it  is  practically 
impossible  to  obtain  any  detailed  information  below  about 
0.2  eV.  It  was  therefore  necessary  to  complete  the  avail¬ 
able  information. 

In  this  study,  optical  constants  based  on  the 
analysis  of  reflectance  measurements  extending  to  225  pm 
and  a  Drude  fit  to  reflectance  extending  to  1000  pm  have 
been  derived.  The  instruments  used  in  the  experiment,  along 
with  the  near  normal  reflectance  measurements  will  be  pre¬ 
sented;  then  the  data  analysis  will  be  discussed. 

Instrumentation  and  Measurements 

Spectrophotometers  to  which  suitable  reflectometers 
(Steyer,  1974)  have  been  attached  were  used  for  the  measure¬ 
ments.  A  Cary  14  was  used  in  the  visible  and  in  the  near 
IR,  a  Perkin  Elmer  398  in  the  IR,  and  a  Beckman  IR  11  for 
the  far  IR  (see  Table  7  )  . 

The  Cary  14  and’ the  Perkin  Elmer  398  were  used  in 
a  double  beam  mode,  and  the  Beckman  11  in  a  single  beam 
mode.  The  Beckman  11  was  modified  to  use  a  liquid  He 
cooled,  doped  germanium  bolometer  as  a  detector.  The  signal 
obtained  was  amplified  by  a  PAR  model  28  lock-in  amplifier 
and  then  sent  to  a  digital  voltmeter  or  chart  recorder 
(Rathmann,  1981). 

The  reflectometers  were  designed  to  give  a  near 
normal  incidence  reflectance.  The  one  used  on  the  Cary  14 


Table  7,  Spectral  region  and  .instrument 


From  Rathm&nn  (1981} 


Region 

Wavelength 

(Frequency) 

Instrument 

Spectral 

Resolution 

Far  IR 

250  ym  -  20  ym 
(40  -  500  cm"*) 

Beckman  IR  11 

%  5  cm 

Mid  IR 

25  ym  -  2.5  ym 
(400  -  4000  cm-1) 

Perkin  Elmer  398 

'v  5  cm 

Near  IR 
to  Visible 

2  ym  «  0. 33  utn 
(5000  -  3.3  x  10K  cm" 

Cary  14 

»> 

MO  Jl 

and  the  Perkin  Elmer  398  is  built  by  Barnes  Engineering 
Company;  the  far  IS  one  was  designed  and  built  in  our  lab. 


An  aluminized  glass  reference  was  used.  Good 
alignment  was  always  checked  before  any  measurement;  the 
aluminum  standard  or  the  sample  was  positioned  and  the 
mirrors  of  the  reflectometer  rotated  until  a  maximum  signal 
was  obtained. 

The  measurement  procedure  was  as  follows : 

1.  Run  the  100%  reflectance  line,  say  R^OO'  witJl 
the  aluminum  standard  in  place. 

2.  Run  the  zero  reflectance  line,  say  RQ,  without 
any  sample. 

3.  Run  the  reflectance  of  a  well-polished  sample, 
say  Rg. 

Prom  the  spectrum  obtained,  it  is  easy  to  find; 


R'tt) 


VR0 

R100~R0 


where  R *  C A )  is  a  relative  reflectance;  this  data  can  be 
used  in  the  Kramers-Kronig  analysis  computer  program  where 
it  is  converted  into  an  absolute  reflectance. 

Reflectance  from  \  ^  0.35  ym  up  to  220  ym  has  been 
measured  and  good  agreement  obtained  with  Arakawa's  pub¬ 
lished  results  in  the  0.'  ym  to  2.5  urn  wavelength  range 
which  the  experiments  have  in  common.  It  was  thus  not 


necessary  to  pursue  measurements  toward  wavelengths  shorter 
than  0.35  ym. 

Data  Analysis  and  Results 

For  use  in  the  Kramers-Kronig  analysis,  reflectance 
was  extended  toward  long  wavelengths  by  curve  fitting.  Free 
electron  behavior  was  assumed  for  glassy  carbon  at  larger 
wavelengths  than  the  last  experimental  point,  and  the  Drude 
model  used  in  the  fitting.  The  parameters  used  were: 
cop  -  2000  cm"1;  ^(1)  »  0,  and  y  a  900  cm"1.  Reflectance 
data  going  from  X  0.02  ym  to  about  106  ym  were  constructed 
by  combining  Williams  and  Arakawa's  (1972)  results,  experi¬ 
mental  results  of  this  study,  and  extrapolated  results,  for 
use  in  the  Kramers-Kronig  analysis  (see  Fig.  21.  ).  The 
results  of  the  analysis  are  presented  in  Figs.  30.  to  32. 
and  in  Appendix  B. 

We  then  proceed  to  use  the  Drude- Loren tz  model  and 
fit  the  reflectance  data  obtained  for  the  Kramers-Kronig 
study;  the  result  is  shown  in  Fig.  33  ,  and  the  parameters 

used  are  presented  in  Table  8. 

12.4  Discussion  of  Results . 

The  reflectance  data  for  graphite  in  the  E  1  C 
orientation  and  that  of  glassy  carbon  have  been  analyzed 
by  the  Kramers-Kronig  method  and  by  fitting  a  Drude-Lorentz 
model  to  the  reflectance  curves.  The  reflectance  data  for 
graphite  in  the  E  //  C  orientation  has  only  been  fitted 


:rt  of  the  dielectric  function  of  glassy  carbon,  as  computed 
Kramers-Krouig  method 


Table  8,  Oscillator  parameters  for  the  fit  of  the  reflec¬ 
tance  curve  of  glassy  carbon 


Parameter 


Oscillator  1 


Oscillator  2 


Oscillator  3 


Oscillator  4 


Value  (cm*1) 


14  , 

.5 

x  10 

5 

X 

103 

9 

X 

10  3 

50 

X 

10  3 

30 

X 

103 

43 

X 

10  3 

120 

X 

103 

110 

X 

10  3 

60 

X 

10  3 

mwMM 


because  the  spectral  range  of  the  data  was  not  large  enough 
to  allow  the  use  of  the  Kramers -Kronig  method. 

In  all  the  fits,  the  oscillator  strength  parameters 
(u£j)  have  been  checked  against  the  sum  rule  (see  Kittel, 
1976,  p.  351  and  Wooten,  1972,  p.  72): 

>* 

we2{w)  dw  ■  l  ui.  -  h  irojp  .  (4,1) 

o  j  3 

The  plasma  frequency  is  given  by  (see  Kittel, 
1976,  p.  289) : 


u 


P 


4trne2 

m 


(42) 


where  n  is  the  density  of  electrons,  e  and  m  the  charge  and 
the  mass  of  an  electron  in  c.g.s.  units.  With  the  assumption 
that  all  the  four  valence  electrons  of  carbon  determine  u>p , 
one  obtains: 

>57r  ajp  «  20.3  x  1012  sec"2 

The  value  of  jj  (u2j  for  each  case  is  givan  in  Table  9. 

In  each  of  the  three  cases 

l  Wpj  <  hir  oj2 
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Table  9.  Values  of  I  . 

j  Pj 


Material 


J  (sec-*) 

'v  4.4  x  10si 
*  3.3  x  10** 

J  2 


Graphite:  £.±  C 
Graphite:  E  ft  C 
Glassy  carbon 


*  6.2  x  10 


This  result  is  expected,  since  the  maximum  possible  value 
for  n  has  been  chosen  for  substitution  in  Eq.  (42)  .  Some 
valence  electrons  might  be  sufficiently  bound  that  the  sum 
rule  (Eq.  41)  is  not  satisfied  over  the  frequency  range  of 
the  experiment. 

The  discrepancy  between  the  experimental  and  fitted 
curves  observed  between  X  ■  0.3  and  15  urn  on  Fig.  25.  is 
due  to  interband  transitions  which  come  into  play  with  in¬ 
creasing  energy.  These  transitions  gradually  increase  the 
reflectance  without  showing  sharp  peaks.  The  Lorentz  model 
used  for  the  fit  cannot  account  well  for  this  phenomenon. 
The  fit  of  the  peak,  at  X^45  ym  resulted  in  the  sharp  drop 
in  the  calculated  reflectance. 

The  band  structure  of  graphite  as  calculated  by 
Mallett  (1981)  shows  narrow  band  gaps  at  X  a  0.08  and  0.25 
ym,  which  account  for  direct  interband  transitions,  respon¬ 
sible  for  the  sharp  peaks  observed  at  those  wavelengths. 
These  peaks  were  easily  fitted. 

Examination  of  Figs.  29  and  33  shows  that  to 
some  extent,  the  reflectance  of  the  E  H  C  polarization  for 
graphite  and  that  of  glassy  carbon  have  the  same  behavior 
as  the  reflectance  of  the  El  C  orientation  for  graphite. 
This  behavior  makes  the  fit  difficult  between  X  ■  0.3  and 
15  ym,  but  the  difference  between  the  experimental  reflec¬ 
tance  and  the  fit  is  unimportant  because  of  the  low 


reflectance  of  the  materials  considered. 
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Results  obtained  by  the  two  methods  of  analysis  for 
glassy  carbon  are  compared  in  Table  10  for  several  wave¬ 
lengths  . 

Comparing  this  Kramers-Kronig  analysis  with  pub¬ 
lished  data,  one  finds  that  for  the  E  1.  C  orientation  of 
graphite,  the  optical  constants  derived  agree  well  with 
those  of  Taft  and  Philipp  (1965)  (see  Table  )  and  Green¬ 
away  et  al.  (1969).  At  the  peak  around  X  ^  0.089  pm,  our 
value  of  9.3  for  e2  is  higher  than  that  of  Taft  and  Phil¬ 
ipp's  (.1965)  value  of  ^  7.0.  The  difference  is  due  to  the 
values  of  the  parameters  used  in  the  short  wavelength 
extrapolations  for  the  Kramers-Kronig  analysis  (see  p.  34) . 

The  optical  constants  obtained  for  the  E  /'  C  orien¬ 
tation  of  graphite  are  to  our  knowledge  the  only  consistent 
ones  available  at  this  writing. 

The  glassy  carbon  optical  constants  agree  well  with 
Williams  and  Arakawa's  (1972)  results  and  extend  them  beyond 
X  10  pm.  Differences  exist  between  our  results  and  those 
published  for  varouis  amporphous  carbons  as  one  can  readily 
see  in  Fig.  is.  One  exception  is  the  result  of  Pluchino 
et  al.  (.1980)  (see  Fig.  32  ).  The  good  agreement:  in  this 
case  is  worth  pointing  out,  even  though  the  comparison  is 
possible  for  one  wavelength  only.  Our  reflectance  measure¬ 
ments  were  performed  on  a  homogeneous  bulk  of  glassy  carbon 
sample,  Pluchino  et  al.  measured  the  scattered  radiation 
by  a  single  particle,  small  enough  to  be  homogeneous.  It  is 


Table  10.  Results  for  index  n  and  k  obtained  by  the 

Kramers-Kronig  method  and  by  oscillator  fit, 
in  the  case  of  glassy  carbon 

The  indices  "OF"  and  MKK"  stand  for  "Oscillator 
Fit"  and  "Kramers-Kronig,"  respectively. 


Wavelength  (ym) 


0.08 

0,15 

0.25 

1.00 

5.00 

10.00 

100.00 

nOF 

1.03 

1.46 

1.40 

1.95 

2.97 

2.88 

5.91 

nKK 

1.04 

1.54 

1.40 

2.15 

3.04 

3.56 

6.09 

kOF 

0.91 

0.38 

0.56 

0.70 

1.15 

1.91 

5.26 

0.95 

0. 19 

0.74 

0.90 

1.03 

1,54 

4.85 

Table  11.  n  and  k  from  this  work  and  those  of  Taft  and 
Phillip  (1965)  and  Philipp  (.1977)  at  the  same 
wavelengths 


Wavelength  (yim) 


VT 

.4 

T5 — 

-  ■ 

1777” 

— 2 - 

5 

This  work 

n  2,5 

2.62 

2.61 

2.7 

3.68 

3.8 

5.2 

16 

k  0.8 

1.28 

1.33 

1.4 

2.4 

2.5 

3.8 

22 

2.6 

3.78 

3.9 

5.1 

15.8 

l 

I 

f 


0.75  1.35  1.36  1.46  2.49  2.6 


4.0 


22.8 


strongly  believed  that  the  agreement  of  our  results  is  due 
to  the  homogeneity  of  the  materials  on  which  the  experiments 
were  performed.  To  find  optical  constants  of  particles, 
many  workers  analyze  reflectance  data  of  particles  pressed 
into  pellets.  Results  obtained  by  this  technique  are  rarely 
in  good  agreement  for  different  samples.  Chylek  et  al. 

(1981)  argued  that  the  differences  are  essentially  due  to 
the  differing  mass  concentration  of  the  particles.  One 
reason  we  have  chosen  graphite  and  glassy  carbon  is  because 
these  materials  are  relatively  homogeneous. 


13.  SMALL  PARTICLE  EXTINCTION 


This  chapter  deals  with  the  extinction  of  small 
carbon  particles.  The  production  of  the  particles  and 
their  characterization  are  presented,  and  the  sample 
preparation  for  the  extinction  measurements  shown.  Then 
tho  measurements  leading  to  the  mass  calibration  and  to 
the  extinction  are  described.  Finally,  results  are  pre-’ 
sented  and  discussed. 

13.1  Particle  Production  and  Characterization. 

Graphite  and  three  other  types  of  carbon  have  been 
investigated.  These  powders  are: 

1.  very  fine  graphite  (Ashbury  #250,  Dixon), 

2  .  standard  lampblack , 

3.  particles  {for  solar  energy  absorption)  produced 
by  Hunt  (1979)  by  acetylene  pyrolysis  at  the 
Lawrence  Berkeley  Laboratory;  those  particles 
are  referred  to  in  this  work  as  "LBL  Smoke", 

4.  particles  were  made  by  evaporating  carbon  in  an 
inert  gas  atmosphere. 

The  chamber  in  which  these  particles  were  made  was 
evacuated  to  2  x  10“ "  torr  with  an  oil  diffusion  pump.  The 
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system  was  then  purged  with  helium  to  reduce  the  partial 
pressure  of  oxygen,  and  reevacuated.  Finally  the  chamber, 
isolated  from  the  pimping  system,  was  filled  with  helium, 
and  an  arc  was  struck  between  carbon  electrodes  (glassy 
carbon) .  The  particles  were  collected  on  appropriate  sub¬ 
strates. 

The  particles  were  observed  with  a  "Scanning  Elec¬ 
tron  Microscope"  (SEM)  and/or  a  "Transmission  Electron  Micro 
scope"  (TEM) .  For  this  purpose  a  small  amount  of  "LBL 
Smoke,"  lampblack,  or  graphite  powder  was  added  to  methane 
and  the  mixture  dispersed  ultrasonically .  The  cloudy  liquid 
was  then  put  in  a  syringe  adapted  to  a  Nueleopore  membrane 
filter  holder  and  filtered.  The  filter  retained  the  indi¬ 
vidual  particles  with  diameter  larger  than  its  pores  and 
also  chain-like  aggregates  of  particles.  Samples  prepared 
with  this  technique  were  observed  with  the  SEM.  Results 
shown  in  Figs.  34 «  and  35*  confirmed  what  is  commongly 
observed:  spherical  shape  for  most  carbon  and  plate-like 
shape  for  graphite  (Walker,  1963).  The  nucleopores  used  had 
0,25  um  diameter  pores.  Many  of  the  particles  presumably 
had  smaller  diameter  than  0.25  um. 

Smoke  made  by  arc  evaporation  of  bulk  glassy  carbon 
was  collected  directly  on  grids  for  observation  with  TEM. 
Results  showed  that  in  this  case,  too,  the  particles  are 
rather  spherical,  although  not  as  perfect  as  observed  for 
other  particles  (see  Fig.  36  ) -  The  heat  generated  by 


15KU  15  5KX 


1U  0023 


Graphite  (Dixon, HPN-2) 


v£  ? 1  <  •  .v ,  4  ’"*< 


::1>. 


r* 


rn 


:r 


Jr1#  '*$Sa 

ta 


JfL1^  »  ,  -A  I 


>  1 


r 


v  \ 


:CV 


4Uf  :*i  W-:S  '*w-*~r  to 

13KU  0? ..  9KX  1U  001  3 


Graphite  (Ashbury,  Micro #250) 


.Min  i  in)  t'li’iMiuii  m  i  r  1 1  n’li’i  i|  " 1  V  i  i'W  of  tw> 
1|  l  .1 1  'll  I  t  I  *  I  ).-|  r  t  i  r  1  cs 


electrons  striking  the  samples  during  the  observation  pro¬ 
cess  is  probably  responsible  for  the  ragged  edges  (Walker, 
1963)  * 


X-ray  powder  patterns  and  electron  diffraction  pat¬ 
terns  observed  in  the  transmission  electron  microscope 
showed  diffuse  bands  with  little  distinction  between  "LBL 
smoke,"  lampblack,  and  glassy  carbon  smokec  In  contrast, 
the  graphite  showed  distinct  diffraction  patterns.  It  was 
concluded  that  the  first  three  are  mainly  amorphous  while 
the  graphite  is  highly  crystalline. 

13.2  Sample  Preparation. 

The  transmission  measurements  from  which  the  extinc¬ 
tion  was  determined  required  suspension  of  the  particle  in 
one  way  or  another.  When  they  were  collected  on  a  usable 
substrate  such  as  quartz  or  LiF,  no  further  sample  prepara¬ 
tion  was  necessary;  otherwise,  special  particle  suspension 
techniques  had  to  be  used. 

For  extinction  measurements  from  the  visible  to  the 
raid  IR,  the  KBr  pellet  technique  was  used.  A  small  amount 
of  the  particulates  (typically  50  to  200  yg)  was  mixed  with 
0.5  mg  of  IR  quality  of  KBr  in  a  vial  that  contained  a  small 
steel  ball.  The  mixture  was  first  hand  shaken  and  then 
shaken  for  5  to  10  minutes  on  a  "Wig-L-Bug"  dental  amalga¬ 
mator.  The  mixture  was  then  transferred  to  an  evacuable 
die  connected  to  a  vacuum-pumping  system  and  placed  on  a 
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hydraulic  press.  After  evacuating  for  about  two  minute's, 
pressure  was  applied  and  slowly  increased  up  to  about 
2000  lbs  per  square  inch.  The  vacuum  was  maintained  for  two 
more  minutes,  then  released,  as  air  was  slowly  admitted  to 
the  system.  Finally,  the  pellet  was  carefully  extracted. 

For  the  far  IR  extinction  measurements,  polyethylene 
powder  was  used  as  the  dispersion  medium;  0.1  mg  of  poly¬ 
ethylene  was  used  in  each  vial.  In  this  case,  shaking  a 
vial  caused  the  mixture  to  stick  c,o  its  wall  because  of 
electrostatic  forces.  To  circumvent  this  problem,  acetone 
was  added  to  the  mixture.  The  vial  was  hand  shaken,  slight¬ 
ly  warmed,  and  left  open  overnight  to  allow  the  acetone  to 
evaporate.  The  mixture  was  then  scraped  off  the  container 
walls  and  transferred  to  the  KBr  die.  The  die  was  heated 
on  a  hot  plate  during  its  evacuation.  A  gentle  hand  pres¬ 
sure  was  exerted  on  the  plunger  of  the  die  at  the  beginning 
of  the  heating.  Still  under  evacuation,  the  die  was  removed 
from  the  hot  plate  and  slowly  cooled.  The  pellet  was  then 
carefully  removed,  with  the  help  of  a  razor  blade,  and 
weighed.  If  m,  is  the  mass  of  the  polyethylene,  m2  the 
mass  of  the  particles  with  which  the  sample  preparation  was 
started,  and  mj  the  mass  of  the  pellet  obtained,  then; 


is  approximately  the  mass  of  particles  in  the  pellet. 
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13.3  Measurements . 

The  ins truments used  to  acquire  data  have  been 
described  earlier  (see  Chapter  12).  With  the  Cary  14,  a 
zero-absorption  trace  was  recorded  with  two  identical  sample 
holders  in  the  two  beams.  The  sample  was  then  inserted  in 
the  front  beam  and  a  reference  blank  was  placed  in  the  back 
beam.  The  spectrophotometer  recorded  the  optical  density 
(OD)  as  a  function  of  the  wavelength.  The  measuring  tech¬ 
nique  on  the  Perkin  Elmer  was  the  same  as  on  the  Cary  14, 
except  that  the  zero  and  the  100%  lines  had  to  be  recorded. 
The  spectra  obtained  gave  intensities  which  permitted  the 
calculation  of  the  transmission. 

The  Beckman  11  was  used  in  a  one-beam  mode.  The 
transmission  through  the  sample  and  the  reference  blank  was 
measured  separately.  The  absorption  of  the  particles  was 
found  from  the  ratio  of  the  two  transmissions.  The  calcula¬ 
tions  of  volume-normalized  extinction  (a)  requires  knowledge 
of  the  mass  density  c  (mass/area)  of  the  particles.  For  the 
pellets  used  in  IR,  the  mass  of  particles  suspnndod  in  a 
matrix  was  large  enough  for  weighing  with  an  analytical 
balance.  In  the  UV,  the  absorption  of  carbon  particles  is 
so  high  that  a  very  small  quantity  of  particles  must  be 
used  to  get  any  transmission.  In  that  case,  it  was  prac¬ 
tically  impossible  to  weight  the  particles;  an  optical 
calibration  method  was  used  to  fina  a. 
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A  sufficient  quantity  of  carbon  smoke  was  collected 

i 

on  glass  slides  and  the  optical  densities  (OD)  measured 
in  the  visible.  Each  slide  was  weighed  with  and  without 
the  smoke.  Knowing  the  mass  m  of  the  particles  on  .a  slide, 
and  the  area  A  of  the  substrate,  a  *  m/A  was  computed. 
Optical  density  was  plotted  as  a  function  of  c  for  several 
samples  of  varying  thickness,  as  shown  in  Pig.  37„  The 
curve  shows  a  leveling  off  of  OD  with  increasing  a,  presum¬ 
ably  due  to  the  increasing  effect  of  multiple  scattering. 

In  order  to  avoid  the  multiple  scattering  region,  the  slope 
of  the  linear  curve  at  small  o  is  used  to  determine  a  cali¬ 
bration  constant  (OD /a)  at  each  wavelength.  This  calibra¬ 
tion  can  then  be  used  to  determine  mass  densities  from 
optical  transmission  measurements  for  samples  that  are  too 
light  for  conventional  weighing  techniques.  By  calibrating 
the  mass  density  in  the  near  infrared  and  the  visible, 
optical  measurements  on  very  light  samples  (i.e.,  5  ug/cm2 ) 
can  be  extended  into  the  highly  absorbing  regions  of  the  UV. 

13.4  Results  and  Discussion . 

In  this  section,  results  of  extinction  measurements 
on  graphite  and  amorphous  carbon  are  presented  and  discussed 
with  respect  both  to  calculations  and  to  other  workers' 
results.  This  is  done  first  for  X  >  0.4  urn  and  then  for 
X  £  0.4  urn. 
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Variation  of  the  optical  density  (OD)  of  glassy  carbon,  as  a  function 
of  the  mass  density  a  at  various  wavelengths 
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In  Fig,  38.  are  plotted  the  extinction  results  from 
graphite  powders  (Ashbury  micro  1250,  Dixon  HPN-2) .  On  the 
same  figure  are  plotted  two  sets  of  calculations,  one  based 
on  a  continuous  distribution  of  spheres  (CS)  theory,  and  the 
other  on  a,  continuous  distribution  of  ellipsoids  (CDS) . 

Although  some  workers  argued  that  bulk  optical 
properties  may  not  be  appropriate  for  micron-sized  particles 
(Pluchino  et  al.,  19  80),.  the  calculations  were  done  using 
optical  constants  derived  in  Chapter  12  with  the  assumption 
that  the  particles  are  large  enough  for  bulk  optical  prop¬ 
erties  to  be  applicable.  A  discussion  of  this  point  can  be 
found  in  the  review  by  Huffman  (1977,  pp.  212-215). 

A  comparison  of  the  two  calculated  results  shows  e 
substantial  difference  in  strength  of  the  extinction  in  the 
infrared  region  from  about  X  'v  1,0  urn  to  100  um.  This  pro¬ 
vides  a  quantitative  illustration  of  how  shape  affects 
extinction  of  carbon  particles  in  the  infrared. 

Visible  and  Infrared  Range 

Most  of  our  studies  on  graphite  were  done  on  the 
Ashbury  powder,  which  was  available  in  large  enough  quanti¬ 
ties  to  permit  the  preparation  of  many  samples.  Therefore, 
the  discussion  of  graphite  experimental  results  will  be 
based  on  the  Ashbury  sample  only. 

Three  sets  of  results  were  obtained  (see  Fig.  38  ) . 
From  X  'v  0.4  um  to  25  um,  extinction  of  particles  in  KBr 
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Experimental  extinction  for  graphite,  compared 
to  sphere  theory  (CS)  and  to  CDE  theory 
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and  extinction  of  particles  in  polyethylene  are  shown? 
above  X  ^  25  pm,  extinction  of  particles  in  polyethylene 
only  is  presented.  The  polyethylene  sample  used  from 
0.4  ym  to  25  ym  was  obtained  from  a  pellet  made  by  , the 
technique  described  in  the  section  on  sample  preparation, 
then  thinned  by  exerting  a  slight  pressure  on  the  pellet 
placed  between  two  microscope  slides  and  heated.  Comparison 
of  the  two  results  in  the  X  ^  0.4  ym  to  25  ym  range  shows 
differences  which  are  due  to  the  suspension  media  used. 

On  the  one  hand,  the  carbon  particles  mixed  well  with  the 
KBr  powder  during  the  vial  shaking  process?  on  the  other, 
the  polyehtylene  and  the  particles  stick  to  the  wall  of  the 
vial  when  shaken,  as  previously  explained. 

The  observed  disagreement  of  theory  and  experiment 
is  most  likely  due  to  the  inappropriateness  of  the  ellip¬ 
soid  theory.  Figure  34.  provides  evidence  that  the  graphite 
particles  are  in  fact  plate-like  particles. 

The  particle  shape  is  probably  responsible  for  the 
violation  of  the  sum  rule  shown  in  Eq.  (40).  Indeed,  in  a 
very  restricted  wavelength  range  (i.e.,  5  ym  to  50  ym) , 
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exceeding  the  maximum  value  which  is  30. 


Although  this  investigation  was  not  planned  as  an 
illustration  of  the  "shape  effect,"  the  comparison  of  exper 
imental  results  to  one  another  and  to  calculations  has 
clearly  shown  the  importance  of  the  shape  effect  in  extinc¬ 
tions. 


Extinction  of  "LBL  smoke"  and  lampblack  was  also 
investigated.  Experimental  results  and  calculations  based 
on  the  same  theories  as  for  graphite  are  plotted  in  Fig. 
39.  Optical  constants  derived  in  Chapter  32  for  glassy 
carbon  have  been  used  in  the  calculation.  Extinction  of 
fresh  "LBL  smoke"  from  \  ^  0.6  ym  to  15  ym  is  pratically 
the  same  as  predicted  by  the  sphere  calculations.  This  is 
not  surprising  since  Fig.  35.  shows  agglomerated  spheres. 
These  chains  and  clusters  were  probably  dispersed  in  the 
shaking  process. 

Pellets  made  of  "LBL  smoke"  were  measured  fresh 


and  also  two  years  later;  they  gave  almost  exactly  the  same 
results  as  new  pellets  made  of  two  year  old  "LBL  smoke." 
Because  of  the  old  pellet  results,  the  increase  in  "LBL 
smoke"  extinction  has  been  assigned  to  a  change  in  the 
optical  constants  of  the  particles,  rather  than  to  the 
shape  effect.  The  nature  of  this  change  is  not  yet  well 
understood;  a  graphitization  hypothesis  has  been  discounted 
because  x-ray  studies  of  the  two-year-old  particles  did  not 
show  any  sharp  lines. 
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aging  effect  observed  and  comparison  to  sphere 
theory  (CS)  and  to  CDE  theory 
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Figure  40.  shows  that  the  lampblack  extinction  is 
higher  than  the  highest  "LBL  smoke"  extinction.  The  lamp¬ 
black  powder  used  might  have  experienced  the  same  aging 
process  as  the  "LBL  smoke"  and  the  difference  in  extinction 
might  be  due  to  a  difference  in  the  ages  of  the  particles. 

Koike,  Hosegawa,  and  Man&be  (1980)  measured  the 
extinction  of  arc  evaporated  graphite  particles  and  parti¬ 
cles  produced  by  burning  benzene  and  xylene  in  air.  The 
measurements  were  performed  from  A  ^  0.21  ym  to  340  urn  on 
samples  prepared  with  the  same  techniques  as  ours.  They 
found  humps  between  A  ^  5  ym  and  15  pm  (Koike  et  al.,  1980). 
Our  data  confirmed  the  presence  of  these  humps  (see  Figs. 

4.5  and  4.6).  Koike  et  al.  also  found  a  peak  near  A  ^  90  urn. 
Neither  the  data  of  Blea  et  al.  (1970)  nor  ours  shows  this 
peak.  A  broad  peak  found  in  our  calculated  result  is  located 
around  A  ^  35  ym  and  is  more  apparent  in  the  sphere  calcula¬ 
tions  than  in  the  CDE  case  (see  Fig.  38  ) , 

Blea  et  al.  (1970)  measured  absorption  coefficients 
of  black  polyethylene  over  a  broad  spectral  range  (A  %  2.5 
ym  to  3,300  urn).  Although  their  scattered  results  were 
assigned  to  differences  in  grain  size,  mass  concentration 
of  carbon  in  polyethylene,  and  to  the  form  of  the  carbon, 
our  extinction  results  from  one  type  of  carbon  to  another 
are  due  to  differences  in  the  shapes  of  the  particles. 

Table  12  shows  experimental  results  of  this  study 
and  those  of  other  workers  at  some  specific  wavelengths. 
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Fig.  40.  Comparison  of  experimental  extinctions  by 
smoke"  and  by  lampblack 
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Ultraviolet  Range 

For  short  wavelengths  the  Rayleigh  approximation, 
on  which  the  sphere  and  CDE  calculations  are  based,  is 
hardly  valid.  A  swi.cn  to  Mie  calculations  is  then  neces¬ 
sary.  Although  the  particles  are  not  individual  spheres  or 
spherical  at  all  (graphite) ,  the  calculations  can  give  a 
rough  basis  for  the  evaluation  of  the  experimental  results. 

Experimental  results  and  calculations  are  shown  in 
Figs.  41.  and  42.  for  Ashbury  graphite  and  glassy  carbon, 
respectively.  It  might  be  surprising  that  the  glassy  carbon 
extinction  peak  is  higher  than  the  graphite  one.  It  has  been 
shown  (Rathmann,  1981;  Foxvog  and  Roessler,  1978)  that  the 
extinction  out  of  the  Rayleigh  limit  range  is  radius  depen¬ 
dent,  and  that  very  fine  particles  have  a  higher  peak.  The 
graphite  used  was  acquired  as  a  finished  powder  with  rather 
large  particle  sizes.  The  glassy  carbon  particles  were  arc- 
evaporated  in  helium  atmosphere,  a  technique  known  for  pro¬ 
ducing  very  fine  particles  (Rathmann,  1981).  Also  the 
particles  were  lightly  deposited  on  the  substrate  to  limit 
their  agglomeration. 

The  graphite  particles  might  obey  some  kind  of  sum 
rule  by  which  the  loss  of  strength  in  extinction  observed 
between  A  «  0.22  vim  and  0.24  urn  is  compensated  by  the  hign 
extinction  in  the  infrared  region. 

The  strength  of  the  experimentally  found  peak  for 
glassy  carbon  is  rather  close  to  the  peak  calculated 
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by  the  Mie  method  for  graphite  (see  Fig.  42  ) .  Although  it 

is  known  that  glassy  carbon  cannot  be  significantly  altered 
by  heat  treatment  (Halpin  and  Jenkins*  1969),  the  strength 
of  the  peak  is  a  good  basis  for  arguing  that  the  particles 
have  graphitized  in  their  production  process. 


14.  SUMMARY  AND  CONCLUSIONS 


Reflectance  data  for  graphite  Cor  E  1  C  orientation 
and  for  glassy  carbon  were  compiled  from  the  literature  and 
augmented  by  our  own  measurements  and/or  extrapolations, 

They  were  analyzed  by  the  use  of  the  Kramers-Kronig  method 
to  yield  consistent  optical  constants  over  a  wide  spectral 
range.  For  glassy  carbon,  our  results  are,  at  this  writing, 
the  only  ones  available  beyond  X  ^  10  um. 

It  was  also  shown  that  the  homogeneity  of  the 
material  on  which  measurements  (reflectance,  extinction, 
etc.)  are  performed  is  an  important  factor  in  obtaining 
dependable  results. 

Optical  constants  obtained  for  graphite  for  the 
E  //  C  orientation  by  optical  and  by  electron  energy-loss 
methods  over  a  broad  (0.05  to  100  um)  spectral  range  have 
been  combined  to  give  more  realistic  results  than  could 
be  provided  by  each  method  separately. 

Fitting  tne  reflectance  data  of  graphite  and 
glassy  carbon  to  a  Drude-Lorentz  model  was  complicated  by 
the  existence  of  interband  transitions  which  increase  the 
reflectance  without  showing  sharp  peaks.  But  except  for 


the  case  of  E.l  C  polarization  of  graphite,  the  difference 
between  experiment  end  fit  was  unimportant  because  of  low 
reflectances;  good  results  were  thus  obtained  with  this 
method  over  the  spectral  range  covered. 

Differences  exist  between  experimental  and  calcu¬ 
lated  extinctions  by  small  carbon  particles i  High  experi¬ 
mentally  observed  extinctions  are  either  due  to  a  shape 
effect  as  observed  with  graphite,  or  to  an  aging  process 
as  noticed  with  the  other  types  of  carbon  particles.  The 
preparation  of  small  particles  by  arc  evaporation  of  bulk 
glassy  carbon  tends  to  graphitize  the  particles,  thus  en¬ 
hancing  the  extinction.  This  affect  was  observed  around 
\  ^  0.22  iim  for  particles  which  were  nominally  glassy 
carbon  (see  Fig.  42  ) . 

The  poor  agreement  between  experiment  and  calcula¬ 
tions  leads  to  the  conclusion  that  current  models  are  not 
adequate  for  predicting  extinctions  by  small  carbon  parti¬ 
cles  over  a  large  spectral  range.  However,  the  empirical 
results  provided  by  this  work  should  be  useful  in  designing 
project?  like  the  solar  energy  collector  conceived  by 
Hunt  (1979). 

More  experimental  work  should  be  done  in  order  to 
reconcile  the  optical  constants  for  graphite  for  E  //  C 
orientation,  as  measured  by  optical  and  by  electron  energy- 
loss  techniques.  The  aging  process  and  the  graphitization 
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of  the  particles  observed  in  this  work  need  to  be  studied 
further  and  quantified. 

Progress  in  reconciling  experimental  and  theoretical 
extinctions  by  small  carbon  particles  will  probably ' depend 
on  how  well  changes  in  the  particles  are  understood  and 
incorporated  in  calculations;  it  will  also  depend  on  the 
existence  of  theories  which  can  properly  deal  with  the  shape 
of  the  particles. 
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APPENDIX  B 


Table  B-l.  Optical  constants  of  graphite  for  the  E  1  C 
polarization 


X  (urn) 

N 

K 

0.050 

0.381 

mm 

0.055 

0.390 

mXBm 

0.060 

0.400 

0.850 

0.065 

0.500 

1.150 

0.070 

0.650 

1.330 

0.075 

0.850 

1.700 

0.080 

1.  000 

2.000 

0.083 

1.900 

2.450 

0.090 

2.600 

1.650 

0.095 

2.630 

1.350 

0.100 

2.400 

0.900 

0.110 

2.  150 

0.500 

0.120 

1.900 

0.130 

0.125 

1.700 

0.100 

0.140 

1.470 

0.130 

0.150 

1.  100 

0.160 

0.160 

0.960 

0.280 

0.180 

0.630 

0.650 

0.200 

0.800 

1.200 

0.220 

0.850 

2.000 

0.240 

1.200 

2.400 

0.250 

1.600 

2.450 

0.260 

2.  100 

2.650 

0.270 

2.550 

2,430 
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1  .700 
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2.620 

1.360 
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2.610 

1*240 

0.450 

2.600 

1.250 
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2.600 

1.370 
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2.700 

1.360 

0.700 

2. 750 

1.480 

o.aoo 

2.850 

1.630 

0.900 

2.92  0 

1.740 

1.000 

1.000 

1.800 
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A  (pm) 


2 
5 

1.600 

2.000 

2.600 

3.000 

3.300 
4.000 

4.300 
5.000 
6.000 
7.000 
8.000 
9.000 

10.000 

12.000 

13.000 

18.000 

20.000 

23.000 

30.000 

35.000 

40.000 

50.000 

60.000 

70,000 

80.000 

90.000 

100.000 


3.700 
3.830 

4.300 
4*300 
4.480 

4.730 
5.030 
3.230 
5.400 
5.400 

5.700 
3.980 
6.100 
6.140 
6.250 
6.550 

6.730 

7.500 

3.500 

10.000 

11.500 

16.500 

20.700 

22.500 

21.500 
20. 000 

13.300 


2.000 
2.180 
2.430 
2.520 
2.750 
2.95  0 
3.330 
3.550 

3.730 
3.830 
4.030 
4.450 
4.950 
3.250 
3.530 
6.650 
G.000 
9.000 

11.000 

13.000 

15.500 

18.300 

20.500 

23.500 

22.000 

19.300 

18.900 

16.800 

21.300 
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Table  B-2.  optical  constants  of  graphite  for  the  E  //  C 
polarization 


X  (pm) 


X  (um) 


0.050 

0.055 

0.060 

0.065 

0.070 

0.075 

O.OSO 

0.085 

0.090 

0.095 

0.100 

0.110 

0.120 

0.125 

o.uo 

0.150 
0.160 
0.180 
0.200 
0.220 
0.2  V0 
0.250 
0.260 
0.270 
0.300 
0.350 
0.400 
0.450 
0.500 
0.600 
0.700 
0.900 
0.900 
1.000 


0.980 

0.950 

0.850 

9.800 
0.790 
0.833 
0.940 
0.960 

1.  30 
1.150 
1.300 
2.200 
2.700 

2.800 
2.550 
2.230 

2.  100 
1.600 
1.220 
1.100 
0.990 
1. 000 
1 .050 
1.120 
1.250 
1.370 
1.400 
1.370 
1.330 
1.300 
1.300 
1.300 
1.300 
1.350 


0.360 
0.420 
0.480 
0.600 
0.780 
0.833 
0.950 
1.05 
1.200 
1.450 
1.600 
1.900 
1.400 
1.000 
0.450 
0*250 
0.190 
0.130 
0.150 
0.200 
0.  120 
0.240 
0.230 
0.100 
0.110 
0,100 
0,100 
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0.100 
0.100 
0.100 
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0.100 
0.100 


1 .200 

1.500 
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3.000 

3.500 
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6.000 
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8  .000 
9.000 

10,000 
12.000 
1 5.000 
18.000 
20.000 
25.000 
30.000 
35.000 
40.000 
50,000 
60 .000 
70.000 
80.000 
90.000 
100.000 


1.400 

1.400 
1.470 
1.500 
1 .600 
1.730 
1*860 
1.960 
2.060 
2.200 
2.350 
2.410 

2.400 
2.398 
2.396 
2.330 
2.250 
2.150 
2.050 
1.930 
1.920 
2.000 
2.080 
2.200 
£.300 
2.430 
2.540 
2.630 
2.700 


0.130 
0.100 
0.100 
0.100 
0.150 
0.150 
0.150 
'  0*150 
0.150 
0*150 
0.150 
0.150 
0.150 
0.150 
0.150 
0.170 
0.200 
0.270 
0.350 
0.380 
0.850 
1.070 
1.240 
1.630 
1.930 
2.15  0 
2.400 
2.550 
2.700 
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Table  B-3.  Optical  constants  and  extinction  coefficients  of  graphite 
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Table  B-4.  Optical  constants  and  extinction  coefficients  of  glassy 
carbon  (continued) 
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Table  B  4.  Optical  constants  and  extinction  coefficients  of  glassy 
carbon  (continued) 
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